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FOREWORD 


This report constitutes the Users Manual that accompanies the 
computer simulation of a space radiator system. It summarizes a two 
year project carried out by the School of Mechanical Engineering at 
the Georgia Institute of Technology, Atlanta, Georgia for the NASA 
Manned Spacecraft Center, Houston, Texas under Contract No. NAS 9-10415. 

The results of the analysis phase of the program are contained in a 
separate report entitled "Space Radiator Simulation - System Analysis". 

A third report submitted under this contract covers the development of 
a simplified system simulation and the initial phases of a system 
optimization procedure. 

The contract was entitled "Study of Design Parameters of 
Space Base and Space Shuttle Heat Rejection Systems". The work was 
monitored by Dr. W. E. Simon of the Power Generation Branch of NASA MSC, 
Houston, Texas, and was carried out by Dr. W. Z. Black and Dr. W. Wulff 
as Co-Investigators and Mr. S. M. Morcos, Mr. S. L. Yao and Mr. R. M. 

Hinson graduate students in the School of Mechanical Engineering under 
the direction of Dr. S. P„ Kezios. Mr. C. R. Caldarale of the Rich 
Electronic Computing Center also contributed to this report. 

The Work carried out by Dr. W. Z. Black is reflected in Section B. 2 
a and b and Section B.3, d through f, of the Permanent Part of the 
program. The program units of the Selective part described in Section 
B.3 were completed by Mr. Morcos under the direction of the Co-Investigators. 

Dr. Black is responsible for the coding of programs No. 14 
through 22, for the supervision of coding programs No. 23 through 30 
which was carried out by Mr. Caldarale and Mr. Hinson and for the super- 
vision of coding programs No. 51 through 60 which was carried out by 
Mr. Morcos. Dr. W. Wulff is responsible for the general program structure 


ii 



described in Section B.l and the program units described in the Sections 
B.3a through c and B.3g of the Permanent Part of the program. Dr. Wulff 
coded programs 1 through 7, 31 through 38 and 40 through 42 and super- 
vised Mr. Yao in the coding of programs No. 8 through 13 and No. 39 
and supervised Mr. Morcos in the coding of programs No. 43 through 50 
and No. 61 and 62. 
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SUMMARY 


This manual describes a computer program that simulates the perfor- 
mance of a space radiator. The program basically consists of a rigorous 
analysis which analyzes a symmetrical fin panel and an approximate analysis 
that predicts system characteristics for cases of non- symmetrical operation. 
The rigorous analysis accounts for both transient and steady state performance 
including aerodynamic and radiant heating of the radiator system. The 
approximate analysis considers only steady state operation with no aerodynamic 
heating. 

The first section of this manual contains a description of the " • 
radiator system and instructions to the user for program operation. The 
input required for the execution of all program options is described. 

Several examples of program output are contained in this sectipn. Sample 
output includes the radiator performance during ascent, reentry and orbit. 

The second section of the user's manual describes the structure 
of the entire program as well as the function of each program unit. All 
interfaces between the various program units are described so that the user 
may remove a single unit or group of program units that have general utility 
for use in other analyses. This section concludes with a program listing. 
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A. USERS INSTRUCTIONS 


1 . Program Objectives and Capabilities 

The program units described in this manual constitute a transient 
heat transfer analysis of a space radiatior heat rejection system. The analysis 
is intended to simulate the radiator system of the space base and the space 
shuttle for conditions of arbitrarily prescribed combinations of aerodynamic 
heating and incident solar, earth albedo and planetary irradiation. 

An exact analysis is performed for the radiator system consisting of 
equally spaced parallel coolant channels existing in one plane connected by 
plane fin panels of trapezoidal cross-section. The midplane between tubes and 
a parallel plane passing through the tube axis are assumed to be planes of 
symmetry for the purpose of the exact analysis. The program results are 
extended over the entire radiator panel on the basis of this symmetry. 

An approximate analysis is used to extend the exact program simulation 
to cover the situation of the non- symmetrical radiator panel. This extension 
predicts on the assumption of one -dimensional analysis, the system performance 
when the radiator panel is curved and when the tubes are unequally spaced, 
non-symmetrically loaded or U-shaped. 

The exact system analysis predicts both transient and steady-state, 
two-dimensional temperature profiles, local and total heat transfer rates and 
coolant fluid temperature and pressure profiles in the flow channel. The 
analysis also predicts the protection layer thickness which covers the coolant 
channel and the inlet and exit manifold tubing. Total system mass and projected 
area are calculated. 

The program units described in this manual serve to simulate the 
radiator system performance. Even though successive runs with system parameter 
variations are feasible, the task of system optimization is dealt with in the 
simplified analysis which is described in detail in the Simplified Analysis 
Manual . 

The purpose of this manual is to describe the program units that make 
up the computer simulation, the necessary input parameters for program execution 
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and the typical program output for simulated radiator operation. The deck 
composition is described in part 2, and the input data set is explained in 
part 3. Part 4 is devoted to a description of the program output and part 5 
gives examples of typical radiator performance for ascent, reentry and orbital 
conditions. Section B describes each subprogram in detail by describing their 
input and output variables as well as the function of each program unit. 

The simulation program is separated into modular units or subprograms 
so that units of general utility may be removed and run separately. This 
arrangement also facilitates program modifications if future changes are 
required. The subprogram concept also reduces the extent of program checking 
when new program units are added to the existing program package. 

The entire program consists of 62 separate program units and the 
Data Set. The source language is FORTRAN V. The program is executed from 
catalogued files, using temporary data files as batch processing on the UNIVA'C 
1108 computer operating under EXEC 8, Version 27.20.143, at the Georgia Tech 
Rich Electronic Computer Center. 

Except for the control statements necessary to collect the program 
into elements to be stored on mass storage for subsequent program execution 
from mass storage, the Executive Control Statements are not discussed. Control 
statements may be freely chosen by the user in accordance with the Executive 
Control Language. All statements are in 026 key punch code. 

a. System Description 

The exact simulation program predicts the performance characteristics 
of the typical coolant channel and fin segment shown in Fig. 1. The fin profile 
has a trapezoidal cross-sectional area. The parallel plane passing through 
the tube axis and the plane bisecting two adjacent tubes are assumed to be 
planes of symmetry. 

The program has property options that allow the user to select from 
a number of different coolant fluids, structural and meteoroid protection 
materials and surface coatings. The materials available for the protection 
layer, tube and fin are beryllium, copper or aluminum in any combination. The 
coolant fluids may be selected from a gas, several liquids or a liquid metal. 
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Fig. 1 Typical Fin Element for Exact Analysis 
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Possible choices are helium, silicon oil, the two f luorochemicals FC-43 and 
FC-75 and the liquid metal NaK. The entire radiator system is covered with 
the same passive thermal coating Z-93 which has optically diffuse but wave- 
length and temperature dependent optical properties. 

The approximate analysis predicts the performance of the non- symmetrical 
radiator panels shown in Fig. 2. This one-dimensional analysis predicts the 
approximate location of the adiabatic plane separating two adjacent coolant 
channels when the radiator panel is curved (Fig. 2a) or the coolant channels 
are shaped in form of the letter U (Fig. 2b). Other situations which the 
approximate analysis is capable of treating involve the radiator panel with 
non-uniformly spaced coolant channels and coolant channels that do not all 
have identical inlet fluid mass flow rates or identical inlet fluid temperatures. 

b. Operating Conditions 

The program simulates radiator characteristics for ascent and reentry 
of the shuttle radiator system as well as orbital conditions for the shuttle 
and space base radiator. Typical system performance for these conditions are 
shown in part 5 of this manual. 

The altitude and velocity profiles of the shuttle vehicle during 
ascent and reentry are supplied by the user. These profiles are read into the 
program from card images. 

The irradiation data to be supplied by the user is subdivided into 
three basic categories: solar, earth albedo and planetary irradiation. This 

data may be either read from card images or from a previously prepared data 
file which contains the irradiation data in the format specified by the MRI 
radiant flux program. 

' * 

The program is capable of system simulation for both transient and 
steady state conditions. Under the steady state option, integration of the 
system differential equations proceeds from the user supplied initial condition 
until steady state is reached. All environmental conditions such as incident 
radiant flux remain constant during the integration process. Under the transi- 
ent option, the environmental conditions are allowed to change with time as 
the integration process proceeds. 




(a) Curved Radiator Panel 



(b) U Shaped Coolant Channels 


Fig. 2 Radiator Panel For Approximate Analysis 


& 


6 


Provision has been made within the program to allow for the variation 
in coolant fluid temperature and mass flow at the inlet plane of the coolant 
channel. These variations must be supplied by the user as sets of discrete 
data values punched on cards. - 
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2. The Deck Composition 

The program structure is dictated by the calling sequence of the 
Runge-Kutta-Simpson integration procedure (see Section III-2 of the System 
Analysis Manual). The Source Deck consists of two parts. The first part, 
referred to as the Permanent Part, requires no attention from the user unless 
he plans major program modifications, for instance to accommodate geometrical 
configurations other than the one considered here. The Permanent Part contains 
the first 42 program units. 

The second part, called the Selective Part, is assembled from 
case to case by the user in accordance with his choice of the 

(i) coolant fluid 

(ii) flow channel material 

(iii) meteoroid protection layer material 

(iv) fin panel material 

(v) thermal control coating 

There are always 20 program units in the Selective Part. 

Following the Source Deck is the Data Set which is discussed in 
Section A-3. It consists of six groups, one for each of the following input 
specifications : 

(i) parameters which specify non-symmetrical radiator panel 

(ii) ascent and reentry profiles 

(iii) solar, albedo and planetary irradiation 

(iv) coolant fluid inlet temperature and mass flux histories 

(v) system specifications 

(vi) program options. 

Most of the user's activities will be reflected in his manipulation of 
the Data Set, especially of the data in groups (v) and (vi) above. 
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3. The Input Data Set 

Input data are divided into six groups and assembled in the order of 
increasing frequency of changes so that data records most likely to be changed 
are at the end of the Data Set. This enables the user to keep large portions 
of the Job Deck on mass storage. The six groups in the Data Set define 

(i) the parameters which specify the non- symmetrical radiator panel 

(ii) the ascent and reentry profiles 

(iii) the incident thermal radiant flux as a function of time 

(iv) the coolant fluid inlet temperature and mass flux histories 

(v) the system specifications, other than material selections 

(vi) the options of program execution. 

The program may be executed successively after reading new data 
records beginning at any of the six data groups listed above. Provisons have 
been also made for restarting the program by reading new data records from one 
group and skipping the following data group or only skipping portions of the 
succeeding data groups. The restart options are determined by the value of 
the restart integer MRSTRT. A description of this integer is given in part 
(f) of this section. 

None of the data groups contains a fixed number of records. The 
number of records in groups (ii), (iii) and (iv) are determined and specified 
to the computer by the user. Groups (i) , (v) and (vi) contain a minimum of 
records as specified below; the number of additional records in these groups 
and of additional sets of groups (ii) , (iii) and (iv) depend on the number of 
additional executions and parameter variations during the same job. In the 
description of data records (card images) the numbering starts from one in 
each group. 
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a) Alphabetical List of Input Data 


All input variables needed for program operation are listed below 
with their units and a brief description of the meaning of the variable. 


VARIABLE 

NAME 

UNITS 

REFERENCED IN 

DESCRIPTION 

ABEMF 


FINA 

ratio of solar absorptivity to 
infrared emissivity of coating 
material 

AL 

ft 

TUBEA 

array of tube lengths 

ALIMIT 

• 

. RUN0PT 

absolute error limit per inte- 
gration step 

ALPHA 

gm/ (day ft 2 ) 

PR0TLR 

constant relating meteoroid flux 
to mass 

ALTA 

ft 

ascent and 
reentry data 

elements in altitude array used 
for ascent altitude profile of 
shuttle 

ALTR 

ft 

ascent and 
reentry data 

elements in altitude array used 
for reentry altitude profile of 
shuttle 

AMAN 

ft 2 

MANIFD 

total manifold area projected 
into plane of fin 

AN 


PROTLR 

see equation 8.6 System Analysis 
Manual 

ATK 


PR0TLR 

see equation 8.8 in System 
Analysis Manual 

BTA 


PR0TLR 

constant relating meteoroid flux 
to mass 

D 

in 

TUBEA 

tube internal diameter 

DITBI 

in 

TUBE 

internal tube diameter 

DTWRTE 

hr 

RUN0PT 

fixed time interval between data 
printout 

GAMMA 


PR0TLR 

constant used to adjust observed 
to predicted meteoroid penetra- 
tion depth 
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H 

in 

FINA 

HFNI 

in 

FIN 

ICASE 


QNML,MRI 

ISYM 


SYSTEM 

ITAPE 


QNMLjMRI 

LFLD 


RUN0PT 

LIMWRT 


RUN0PT 

LTS 


RUN0PT 

LTT 


RUN0PT 

M 

lbm/hr 

FLUIDA 

MD0TI 

lbm/hr 

FL0W 

MRSTRT 


GINPT 

MST0TR 


RUN0PT 

MZ 


TUBE 

NA 


ascent and 
reentry data 

NC0NV 


RUN0PT 


array of fin heights 

fin height from root to tip 

integer value which specifies 
set of irradiation values on 
MRI file 

integer which specifies type 6f 
non-symmetrical loading of- tubes 

integer value to specify card 
input (=0) or MRI file input 
0) for irradiation data 

control integer used to select 
between variable or fixed inlet 
fluid properties 

maximum number of data record- 
ings during integration toward 
steady state 

control integer used to select 
initial system temperature 

control integer used to select 
between temperature dependent or 
independent surface coating 
properties 

array of mass flow rates at inlets 
to tubes 

total coolant mass flow rate 

control integer used for selection 
of reentry point for restarting 
program 

control integer for selecting 
transient or steady state 
computation 

number of nodal points along tube 
axis (<_ 10) 

number of instances for specified 
ascent velocity and altitude data 

control integer for selection of 
orbit, ascent or reentry 
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NFLD 


inlet coolant 
fluid history 

number of mass flow rate and 
temperature data values 

NR 


ascent and 
reentry data 

number of instances for specified 
reentry velocity and altitude 
data 

NRMP 


PR0TLR 

number of radial nodal points in 
the protection layer (<. 5) 

NRTBI 


TUBE 

number of radial nodal points in 
tube wall (<. 5) 

NSRD 


SYSTEM 

number of sides of the radiator 
panel which radiate (1 or 2) 

NT 


SYSTEM 

number of flat symmetrically 
loaded panels 

NTBS 


TUBE 

number of flow channels 

NTM 


QNML ,MRI 

number of instances for specified 
irradiation data (< 100) 

NX 


FIN 

number of nodal points on fin 
perpendicular to flow channel 
(< 10) 

PHI 


PROTLR 

see equation 8.1 System Analysis 
Manual 

PHIN 

degrees 

MRI 

array of angles that relate fin 
panel orientation with MRI data 

PHN 

degrees 

QNML 

angle of panel relative to 
MRI data 

PIN 

lbf/ft 2 

FLUIDA 

array of pressures at inlet to 
tubes 

PO 

lbf/ft 2 

FL0W 

entrance pressure into flow 
channel 

PR0B 


PR0TLR 

probability of no damage caused 
by meteoroid impact 

RH0FNI 

lbm/ ft 3 

FIN 

fin material density 

RH0MET 

g/cm 3 

PR0TLR 

meteoroid density 

RH0MPI 

lbm/ ft 3 

PR0TLR 

protection layer density 
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RH0TBI 

lbm/ ft 3 

TUBE 

tube wall density 

RLIMIT 


RUN0PT 

relative error per integration 
step 

SR00TI 

in 

FIN 

fin root thickness 

STAGX 

ft 

FIN 

distance from stagnation point 
on shuttle to center of radiator 
panel along a streamline " 

STBI 

in 

TUBE 

tube wall thickness 

STIPI 

in 

FIN 

fin tip thickness 

TA 

sec 

ascent and 
reentry data 

elements in time array used for ' 
ascent velocity and elevation 
profiles 

TAU 

days 

PR0TLR 

time radiator is exposed to 
meteoroid environment' 

TEND 

hr 

RUN0PT 

termination time for transient 
performance calculation 

TH 

in 

FINA 

array of fin thickness 

THETA 


PR0TLR 

see equation 8.1 System Analysis 
Manual 

TI 

R 

RUN0PT 

initial system temperature 

TIFLD 

R 

coolant fluid 
inlet temper- 
ature history 

array of coolant fluid inlet 
temperatures 

TIN 

R 

FLUIDA 

array of temperatures at inlets 
to tubes 

TO 

R 

FL0W 

entrance temperature into flow 
channel 

T0 

min 

QNML,MRI 

time in minutes at which 
irradiation data is required 

TR 

sec 

ascent and 
reentry data 

elements in time array used for 
reentry velocity and elevation 
profiles 

TSTAR 

R 

SYSTEM 

array of sink temperatures 
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VELA 

ft/sec 

ascent and 
reentry data 

elements in velocity array used 
for ascent velocity profile of 
shuttle 

VELM 

ft/sec 

PR0TLR 

average meteoroid velocity 

VELR 

ft/sec 

as cent and 
reentry data 

elements in velocity array used 
for reentry velocity profile of 
shuttle 

VERTX 

ft 

FIN 

total dimension of radiator 
panel in direction parallel to 
acceleration of gravity 

W1 

in 

TUBEA 

array of distances between side 
segments of tubes when U-shaped 

W2 

in 

TUBEA 

array of distances between base 
segments of tubes when U-shaped 

WIFLD 

lbm/hr 

coolant fluid 
inlet mass 
flow rate 
history 

array of coolant fluid inlet mass 
flow rates 

XL 

ft 

TUBE 

tube length 
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b) Parameters Specifying Non- Symmetrical Radiator Panels 

Normal program operation assumes that radiator conditions are 
symmetrical about each coolant fluid channel. Several factors may destroy 
this symmetry. One may be non-uniform spacing of the tubes. Another factor 
may be different inlet temperatures or mass flow rates into adjacent tubes. 

A third factor could be attributed to a curved panel resulting in a variation 
to sink temperature between adjacent fin segments. These factors are speci- 
fied in the NAMELIST' s that follow. 

Numerical values for a non-symmetrical panel, for irradiation data 
and for system specifications are read into the computer in NAMELIST format s 
one NAMELIST each with name and input data as follows: 

i) for non-symmetrical panel specifications (see Section (b)) 

SYSTEM for selection of system parameters when panels and/or 
flow conditions are non-symmetrical 

TUBEA for non-symmetrical tube specifications 

FINA for non-symmetrical fin specifications 

FLUIDA for non-symmetrical fluid specifications 

MRI for non-symmetrical irradiation data 

ii) for irradiation specifications (see Section (c)) 

QNML for irradiation data under symmetrical radiator panel 
conditions 

iii) for system specifications (see Section (f)) 

TUBE for flow channel specifications 
FL0W for coolant fluid flow specifications 
FIN for fin parameter specifications 
PR0TLR for the protection layer specifications 
MANIFD for the manifold input data. 

Each input data record under NAMELIST format has 

(i) a blank in its first column, 

(ii) a $ - sign in its second column, 

(iii) the NAMELIST name, starting in the third column, followed by 
one or more blanks. 
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(iv) the string of data in free field format as follows: 

A = b , 

where (A) is the variable name and (b) is a number; imbedded 
blanks between (A), (=) , (b) , ( , ) and a following (A) are 
permissible, 

(v) the comma of the last datum is to be placed by a $ - sign 
which delimits the NAMELIST data string. 

For more details concerning NAMELIST formating, the reader may consult standard 
texts on FORTRAN PROGRAMMING. The NAMELIST data strings must be sequenced in 
the order in which they are listed above. 

i) NAMELIST/ SYSTEM / contains 4 variables 

NT number of flat symmetrically loaded panels, each 

panel containing NTBS symmetrically loaded tubes 
(see Fig. 3 ) 

NSRD = 1 for one side of the panel radiating 

= 2 for both sides of the panel radiating 

Note: if panel is curved NSRD is set equal to 1 

ISYM = 0 for symmetrical panel in one plane and straight 
parallel equally spaced tubes 

= 1 for non-symmetrical curved panel consisting of 
straight fin segments and straight parallel tubes 

= 2 for non-symmetrical panel in one plane with U- 
shaped tubes 

= 3 for non-symmetrical curved panel consisting of 
straight fin segments and U-shaped tubes 

= 4 for non-symmetrical panel in one plane with 
straight parallel tubes but non-symmetrical 
loading 

TSTAR an (NT, 3) array of sink temperatures in °R. 

These values are ignored if MRI file is used 
for incident flux data (see Section (d)). 

Also: It' the radiator panel is flat only (1, NSRD) 
of the array is required. 

The following four NAMELIST' s are not read and should not be placed in the 
program if ISYM =0. If ISYM ^ 0 these four NAMELIST's should appear in the 
order shown. 
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ii) NAMELIST/TUBEA / contains 4 variables 

D the single value of all tube internal diameters 

in inches 

AL an array of NT tube lengths in feet. These 

values are used only for straight tubes. For 
curved tubes, the lengths are calculated from 
the values of W1 and W2 which are read in below 

W1 distance between U-shaped tubes along inlet and 

exit portion of the tube in inches (i.e., both 
side segments of the U) NT + 1 values 

W2 distance between U-shaped tubes along bottom 

segment of U in inches, NT + 1 values (see Fig. 

4 for details of dimensions W1 and W2) 

iii) NAMELIST/FINA / con tains 3 variables 

H NT + 1 values of the fin height in inches from 

tube centerline to midplane between tubes. This 
value is used on for the case of straight tubes. 
Fin height for cases of the U tube is calculated 
from W1 and W2 in NAMELIST TUBEA 

TH NT values for the fin thickness in inches 

ABEMF ratio of the solar absorptivity to infrared 

emissivity of the surface of the radiator system 


iv) NAMELIST/FLUIDA / contains 3 variables 


TIN NT values for the inlet fluid temperatures in R 

PIN NT values for the inlet fluid pressure in lbf/ft 2 

M NT values for the inlet fluid mass flow rate in 

lbm/hr 


v) NAMELIST/MRI/ contains 5 variables 


ITAPE =0 to read irradiation data from cards 

>0 to read irradiation data from MRI file. Value of 
ITAPE designates the input unit used to read file 
data 

ICASE integer used to designate the case number of 

stored data to be used for input. ICASE = 3, 
for example, would designate the third set of 
data stored on the MRI file to be used for 
current input. If ITAPE =0, this integer is 
ignored. 
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Adiabatic Surfaces 



Fig. 4 The Dimensions W1(I) and W2(I) 
for U Shaped Coolant Channels 
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NTM number of instances at which incident flux values 

on to be specified if either card or file data is 
used; NTM <_ 100. If ITAPE > 0, NTM must be the 
number of data values stored on the MRI file 

T0 time in minutes during orbit at which irradiation 

data is required. Time is measured from zero time 
specified on MRI output data. If ITAPE = 0 this 
value is ignored 

PHIN array of NT values of the angles in degrees 

between the normals to the NT flat radiator 
segments and the reference line used as the zero 
angle reference on the MRI program 

c) Ascent and Reentry Data 

Velocity and altitude profiles for the shuttle vehicle during both 
ascent and reentry are read from cards which have the format listed in the 
table below. Selection of ascent or reentry phases is determined by the 
control integer NC0NV (see Section (e) which explains the Execution Control 
Options) . 


Card No. 

F0RMAT 

Variable 

1 

2110 

NA,NR the number of ordered pairs of 
elements in the velocity-time and 
altitude- time arrays for ascent and 
reentry profiles, respectively 

^ though 
INT (NA/10 + 0.9) 
+ 1 = N 

10F8.3 

TA, elements of time array (NA values) 
in seconds, selected for ascent velocity 
and elevation profiles 

N + 1 
through 

N + INT (NR/ 1040.9) 
= n 2 

10F8.3 

TR, elements of time array (NR values) 
in seconds, selected for reentry velocity 
and elevation profiles 


INT(x) = Integer part of x 
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N + 1 
through 

N 2 P INT (NA/ 8+0 . 875) 

= n 3 

8E10.3 

VELA, elements of velocity array (NA 
values) in ft/sec, selected for ascent 
velocity profile of orbiter 

Njf 1 
through 

N 3 + INT (NR/ 840 . 8 75) 
N 4 

8E10.3 

VELR, elements of velocity array (NR 
values) in ft/sec, selected for reentry 
velocity profile of orbiter 

N iff- 1 
through 

N 4 4- INT (NA/ 840. 875) 
= n 5 

8E10 . 3 

ALTA, elements of altitude array (NA 
values) in ft, selected for ascent 
elevation profile of orbiter 

n 5 + 1 

through 

N 5 4- INT (NR/ 840. 8 75) 

8E10.3 

ALTR, elements of altitude array (NR 
values) in ft, selected for reentry 
elevation profile of orbiter 


d) Irradiation Data 

Incident radiant flux on the space base during orbit may be either 
read from card images or it may be read from a user supplied data file 
(MRI program) . The selection of card or file input is determined by the 
value of the control integer ITAPE. If ITAPE = 0, the irradiation data is 
read from cards in the format shown below. If ITAPE > 0, the irradiation 
data is read from unit ITAPE. The data is assumed to be in the format 
specified in the MRI manual "Orbiting Satellite Surface Temperature Predic- 
tion and Analysis," MRI Project No. NAS9-1059. 

The input variables used to input the irradiation data are to be 
supplied in the NAMELIST format called QNML , if ISYM = 0. The input strings , 

must be in the sequence shown below. If ISYM 0, these values are supplied from 
the NAMELIST called MRI 

NAMELIST/ QNML / contains 5 variables 

ITAPE =0 to read irradiation data from cards 

>0 to read irradiation data from MRI tape. Value 
of ITAPE designates the input unit used to read 
file data 

ICASE integer used to designate the case number of 

stored data to be used for input. ICASE = 3, 
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for example, would designate the third set of data 
stored on the MRI file as the irradiation' data to 
be used for current input. If ITAPE = 0, this 
integer is ignored 

NTM number of instances at which incident fluxes are 

to be specified if either card or file data is 
used; NTM < 100. If ITAPE > 0, NTM must be the 
number of data values stored on the MRI file 

T0 time in minutes during orbit at which irradiation 

data is required. Time is measured from zero time 
specified on MRI output data. If ITAPE = 0 this 
value is ignored 

PHN single angle in degrees between normal to panel 

and reference line used as a zero angle reference 
in MRI program. 

If the irradiation data is to be read from cards, it must be in the format 
shown below with the cards immediately following the NAMELIST input cards for 
QNML. 


Card No. F0RMAT Variable 

1 through 10F8.3 On each card; 

NTM (i) time TM in seconds 

(ii) 3 values of solar 

(iii) 3 values of albedo 

(iv) 3 values of planetary 


irradiation 
(Btu/hr ft 2 ) 


Incident at TM and each 
coming from one of three 
directions in this order: 


(i) -30° 
(ii) 0° 
(iii) +30° 


from the panel 
normal (upper side) 


NTM + 1 
through 
2 *NTM 


10F8.3 


Same data as above, now for lower fin 
panel side, only if NSRAD = 2 
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e) Cool^t Fluid Temperature and Mass Flux Histories 

When the temperature and mass flow rates entering the coolant 
channels vary with respect to time, the variation in these quantities are 
entered from cards in the format indicated below. 


Card No. 


F0RMAT 


Variable 


1 

2 through 
NFLD + 1 


110 

3F20.6 


NFLDTA the number of data values of 
temperature and mass flow rate to be 
entered on each card: 

i) TMEFLD the time in hours cor- 
responding to the given tempera- 
ture and mass flow rate 

ii) TIFLD coolant fluid .temperature 
in R at time NFLDTA 

iii) WIFLD coolant mass flow rate in 
lbm/hr at . time NFLDTA 


f) System Specifications 

The system specifications define the radiator system and follow 
immediately after the ascent and coolant fluid temperature and mass flux histories. 

The first four cards in the group of the system specification 
data contain the names of 

(i) the tube material 

(ii) the fin material 

(iii) the coolant fluid ... 

(iv) the meteoroid protection layer material 

in alphanumeric form, beginning with an alphabetic character. Each name is 
punched on a separate card, beginning in the first column and extending no 
further than through the first twelve columns. These names serve to identify 
the printout of results. The user must supply these names in agreement with 
the composition of the Selective Part described in Chapter I. A. The subsequent 
data records contain numerical data. 
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(i) NAMELIST/TUBE/contains 7 variables 


DITBI 

internal (hydraulic) tube diameter, in 

inch , 

STBI 

tube wall thickness, in inch. 


XL 

tube length, in ft. 


RH0TBI 

tube wall density, in lbm/ft 3 , 


MZ 

number of nodal points along the tube 
(no. of intervals plus one), <_10, 

axis 

NRTBI 

number of nodal points in radial direction, in 
the tube wall, <. 5, 

NTBS 

number of flow channels. 



(ii) NAMELIST/ FLOW / con tains 3 variables 

MD0TI total coolant mass flow rate, in lbm/hr, 

TO entrance (and reference) temperature in degrees 

R, 

PO entrance (and reference) pressure in lbf/ft 2 . 

(iii) NAMELIST/FIN /contains 7 variables 

number of nodal points on the fin, perpendicular 
to the channel axis , NX <_ 10 , 

fin root thickness, in inch, 

fin height, from root to tip, in inch, 

fin tip thickness, in inch, 

fin material density, in lbm/ft 3 , 

distance from stagnation point on shuttle vehicle 
to the center of radiator system measured along 
a streamline, in ft, 

total dimension of radiator panel in the direction 
parallel to the acceleration of gravity, in ft. 

(iv) NAMELIST/PR0TLR/ contains 13 variables 


NRMP 

number of radial nodal points in the protection 
layer, <_ 5 , 

RH0MPI 

protection layer material density, in lbm/ft 3 . 

RH0MET 

meteoroid density in g/cm 3 . 

VELM 

meteoroid velocity in ft/sec. 

TAU 

time that vulnerable area is exposed to meteoroid 
environment , in days , 

PR0B 

probability of no damage caused by meteoroid 
impacts, dimensionless, 


NX 

SR00TI 

HFNI 

STIPI 

RH0FNI 

STAGX 

VERTX 
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ALPHA experimental constant that relates meteoroid 

flux and mass (see Eq. 8.4) in gm/(day ft 2 ), 

BTA experimental constant that relates meteoroid 

flux and mass (see Eq. 8.4) dimensionless, 

GAMMA empirical constant used to adjust predicted 

penetration depths to one observed experimentally 
(see Eq. 8.1) dimensionless, 

PHI empirical constant (see Eq. 8.1) dimensionless, . 

THETA empirical constant (see Eq. 8.1) dimensionless, 

ATK empirical constant used to account for spalling 

on a target of finite thickness (see Eq. 8.8) 
dimensionless , 

AN experimental constant that describes penetration 

depth as a function of angle of incidence (see 
Eq. 8.6) dimensionless. 

(v) NAMELIST/MANIFD / contains 

AMAN total manifold area, projected into the fin plane, 

in ft 2 . 

Any one of these variables may be changed for successive program executions 
by entering them in NAMELI S T/ GINPT/ as discussed below. 


g) Execution Control Options 

The user has the choice of options for the execution of the first as 
well as successive computer runs. 

Two data records are considered to be execution controls: RUN0PT 

and GINPT. 


(vi) NAMELIST/RUN0PT/ con tains llvariables 


MST0TR = 1 
= 2 


to compute steady state conditions 
to simulate transient system performance. 


DTWRTE 

TEND 

ALIMIT 


fixed time interval, in hr, between data 
printout during integration, 

termination time, in hr, for transient 
performance calculation, 

absolute error limit per integration step, see 
Eq. 15.3, in Systems Analysis Manual 
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KLIMIT relative error limit per integration step, see 

Eq. 15 . 3 , in Systems Analysis Manual. , 

TI initial temperature, in degrees R, 

LIMWRT maximum number of data recording during integra- 

tion toward steady state, exclusive of initial 
conditions record and steady state record, 

NC0NV =0 no aerodynamic heating 

= 1 ascent 

= 2 reentry, 

LTT = 0 optical properties of the surface coating are 
independent of temperature 

= 1 optical properties of the surface coating are 
dependent on temperature 

= 2 irradiation on fin on tube is zero but subroutine 
QINCID is called, 

LFLD = 1 coolant fluid inlet temperature and mass flow 
rate are constant 

= 2 temperature and mass flow rate of inlet coolant 
fluid is variable, 

LTS = 0 uniform initial system temperature equal to TI 

= 1 initial system temperature taken to be equal to 
previously obtained values. 

(vii) NAMELIST/ GINPT / contains every variable name listed in items 

(i) through (vi) in addition to NSRD and the following control 
integer. 

MRSTRT = 1 new velocity and altitude profiles 

= 2 new irradiation history 

=3 both new velocity, altitude profiles plus new 
irradiation history 

= 4 new coolant fluid inlet conditions 

= 5 new coolant fluid inlet conditions, plus new 
velocity and altitude profiles 

= 6 new coolant fluid inlet conditions, plus new 
irradiation history 

= 7 new coolant fluid inlet conditions, plus new 
velocity and altitude profiles, plus new 
irradiation history 

= 8 the only new input variables are those specified 
in GINPT. 
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There may be as many NAMELIST/GINPT/ records as computing time permits. If 
GINPT is omitted, the program will terminate. 
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4. Results 

a. Information Produced 

The results produced by the program can be divided into three 
categories each of which is listed on a separate output page. All input 
variables that specify the system and select the program options are printed 
on the first page of output (see Fig. 5 ). The second page contains all 
initial line conditions and system parameters which remain constant during 
the integration process (see. Fig. 6 ). The succeeding pages contain the 
current system description printed during the integration process (see Fig. 7 )• 
If the system is non-symmetrical an additional page of output is printed. 

This page (see Fig. 8 ) lists system parameters when the radiator panel is 
curved, the flow channels are not identically loaded, or the flow channels 
are U-shaped. 

b. Output Format 

The first output record consists of the list of all system speci- 
fications and execution control data as listed in Sections 3(d) and 3(e). 

This output appears with all of the variables in the NAMELISTS TUBE, FL0W, 

FIN, PR0TLR , MANIFD and RUN0PT listed in the NAMELIST format. A sample of 
this output is shown in Fig. 5 . 

On the second page of output the initial coolant fluid flow 
properties are printed in non-dimensional form (see Fig. 6 ). This information 
includes non-dimensional pressure, velocity, temperature, and tube wall 
temperature. The reference quantities are listed below the initial line 
conditions. They are inlet fluid pressure, velocity and temperature. Next 
are printed several dimensionless parameters such as the Reynold's number, 
Prandtl number and Nusselt number for the initial coolant fluid flow conditions. 
The quantity DELTA is the ratio of the internal tube diameter to the tube 
length. The relative pressure is given by equation 3.31 in the System 
Analysis Manual. The Biot number of the tube is also printed. 
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Current System Description 
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A listing of system parameters follows the initial line condition. 

The system parameters include the pertinent geometry parameters, the materials 
and the mass of the tube, fin, coolant fluid and protection layer materials. 

In addition, the thickness of the meteoroid protection layer, the total 
system mass excluding the mass of the manifolds and the total projected areas 
are listed at the bottom of the page. 

Following the preliminary results is a series of current system 
descriptions printed out during the integration process (see Fig. 7 )• The 
first items printed are the current real time and the relative time which 
is the ratio of the current time to the time required for a fluid particle to 
pass through the flow channel. If the calculations are for the space base in 
orbit, the words "in orbit" are printed at the top of the page. If the 
calculations are for either ascent or reentry of the space shuttle, the 
current shuttle velocity and altitude are printed along with the current 
atmospheric temperature. For either the sub-orbital or orbital case, the 
cummulative number of integration steps from the initial conditions is 
printed at the top of the page. 

The next data printed are a summary of the heat flux terms and a listing 
of the non-dimensional fin temperature distribution. The reference tempera- 
ture is the inlet fluid temperature and the reference radiant heat flux is 
the black body emissive power at the inlet coolant fluid temperature evaluated 
for a unit length of fin and tube segment. The total radiant rejection and 
the conduction into the fin from both the inlet and exit manifolds for the 
entire radiator system are also printed. The external heat fluxes incident 
on the system that are printed are the infrared irradiation from earth 
albedo, solar irradiation and the convective flux on the radiator system of 
the shuttle vehicle during ascent and reentry. The energy storage rate for 
the entire system is also listed along with the external heat fluxes. 

The final table listed on the page contains coolant fluid properties 
as functions of axial distance as well as the local cooling rates. All 
properties are tabulated in non-dimensional form and the table is preceded 
by a list of the necessary reference quantities. Properties of the coolant 
fluid that are printed are pressure, velocity and bulk temperature. The 
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non-dimensional interior tube wall temperature and the normalized temperature 
of the meteoroid protection layer temperature are also listed as functions 
of axial distance down the tube. The enthalpy rejection of the fluid and 
fraction of the total enthalpy drop are listed as functions down the tube. 

The latter parameter is helpful in determining the effectiveness of the 
various segments of the tube length in their capacity to reject heat. Below 
the reference properties used for the fluid properties are the inlet and 
exit fluid stagnation enthalpies as well as the total change in stagnation 
enthalpy. These quantities are labeled respectively, HI, HE and DH. Also 
for comparison is listed the enthalpy flux at the reference state (T 0 ,P 0 ) and 
the reference velocity W Q . 

The last table of current output data constitutes either the 
system conditions at the specified termination time, TEND in hours, or the 
steady state condition depending on whether MST0TR = 1 or MST0TR = 2, 
respectively. For conditions of steady state, the words "STEADY STATE IS 
REACHED" are printed after the final data printout. 

When input parameters lead to a condition for which the tube 
centerline or the midplane separating two tubes are not lines of symmetry, 
then a single page of output is printed prior to all other output. This 
page contains a written description of tube arrangement on the radiator and 
a printout describing the panel as either flat or curved. The number of 
non- symmetrical tubes and the number of sides of the panel radiating are 
also printed at the top of the page. 

The next output consists of a description of the sink temperature 
data including orientation of the fin segments relative to the reference 
vector used in the MRI program. 

The tube specifications follow the sink temperature data. The 
data listed are the internal tube diameter and the tube lengths when the 
tubes are straight. If the tubes are U-shaped, the spacings at the inlet 
and outlet segments as well as the crossover segments of the tubes are 
listed instead of the tube lengths. 
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The fin specifications follow the tube data. The data listed are 
fin thicknesses and the distance between the tube centerline and the midplane 
separating two tubes. The effective sink temperature for each fin segment is 
also printed under fin specifications. When both sides of the panel radiate, 
two lines of sink temperatures are printed corresponding to the sink tempera- 
ture both radiator sides. 

The coolant fluid specifications follow the fin data. Inlet fluid 
temperature, pressure and mass flow rates are printed for each non- symmetrical 
tube . 

The final data to be printed on the page is a summary of the system 
conditions for all of the non-symmetrical fin segments. The data printed 
includes the fin bulk temperature, position of the adiabatic plane, the fin 
effectiveness, fin tip temperature, fluid outlet temperature and total heat ' 
rejected by the fluid in each of the tubes. 

This single page of output is then followed by the usual system 
printout as described above, repeated once for each non-symmetrical tube. 

c. Output Listing 

Figures 5, 6 and 7 in this section are typical output listings 
produced by the simulation program. Figure 5 provides a listing of all 
input parameters that appear in the NAMELIST format. Figure 6 shows initial 
line conditions and system parameters and Figure 7 shows a typical printout 
of the current radiator system parameters. These three figures are typical 
of program output when the radiator panel is symmetrical about each flow 
channel. If non-symmetrical conditions exist, output similar to that shown 
in Figure 8 is produced. 
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5. Typical Radiator System Performance 

This section summarizes typical program results for expected 
radiator input conditions during ascent and reentry for the shuttle vehicle 
and during both transient and steady state performance of the radiator system 
in a simulated orbit. 

a. Ascent j 

(i) Mission Description j 

The ascent altitude and velocity profiles used for the system 

output shown in Figures 10 and 11 is given in Figure 9. Other input 

parameters were selected as reasonable conditions expected for a 

typical space shuttle radiator system. The total fin area of 124.5 ft^ 

is subdivided by 10 equally spaced 12-foot long flow channels. The total 

fin height between flow channels is 1 ft. The inlet coolant fluid 

2 

temperature and pressure are 200 F and 60 1b^/in , respectively. The 
system materials are aluminum tubes, fins and meteoroid protection 
material with silicon oil as the coolant fluid. Other system 
parameters are printed as input listing which precedes the current 
system parameters shown in Figure 10a. The panel radiates only from 
one side. Incident radiant fluxes are constant during ascent at 
100 Btu/hr ft^ for solar irradiation, 2.0 Btu/hr ft2 for earth albedo. 
There is no infrared earth irradiation. 

(ii) System Response 

The response of the radiator system during ascent is summarized 
in Figures 10 and 11. A selected number of output listings of current 
system parameters during ascent are shown in Figure 10a. Figure 10b 
is a plot of the total heat rejection from the entire radiator surface 
during ascent as well as the outlet coolant fluid temperature. Figure 11 
is a plot of the aerodynamic heating at various times during the ascent 
phase of the mission. 
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Figure 10(a) Continued 
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Figure 10(a) Continued 



f 


l. 






x 






o 


















CM 




r- h- 


0 




X Ll 


it 










GJ <3 






in in 






X X 






X X 

X X 

in 





X X 

1-4 




x w X 





3 3 3 3 

UJ 




t- »- 

r- k 

x 




CD CO 

CD CD 

0 2 




IO -ft 

in 0 

CL < 




o o 

o o 

3 X 




♦ ♦ 

♦ ♦ 

r- UJ 




O O r- io 

in cl 




O 0,0 IO 

X 




o o 

d- ro 

• UJ 




-ft CM 

CM *H 

x >- 




• • 

• • 

o 






UJ * 






H» X 




11 II 

It It 

z »- 






►ft < 










x a > o 

m 




3 UJ Z X 

in 




o x 

3 O 





in ►-♦ 

O »- 





o a 

<3 in 











X 

X 

o 




3 

UJ 

bl 




* 3 

X *■ 

i/I 




X u. 

O UJ 

X 




3 

a i- 

r~ 




-1 Q 

< 

x 




X ill 0 X 

<0 




X 

z 





X ‘jj 

o in 




< Xl>- o 

• • 




JU.<< 

0 O' 




O Z UJX 

-ft O' 




in *-i!X O 

H 





h; 

vO 




r- i- 

• in 





Z ZiZ 





UJ UJ >• >- 



z * 


O 0:0 O 



o * 


*-4 *— 4 O X 



►ft * 


O OiX UJ 



r— * 


z z 

-d Z 



o * 


M H 

< UJ 

in 


UJ * 




M 


"5 * 






UJ ft 




UJ 


x ft 




X 


* 




ti in 


H * 




H* *-• 


< ft 






Uj * 




UJ >• 


X * 

4b 



> H* 


ft 




»-i *-« 


»- * 


>- 


r- o 


z ft 


Ll 


< o 


< ft 




_j _j 


« 


X: 

UJ UJ 


o * 


XX X 

X > 


< ft 



X X 



x * 


XX X 



* 


33 3 



c. # 


H»- *“ 



z ft 

a 

cn:o CD 



< ft 






* 

o 

^ in ro 



z * 

r- 

o 

o o 



o * 

0 

b 

♦ b 

1 

►ft # 

• 

0 co a 

tr H* 


►— ft 

O' 

^ in co 

X X 


3 * 

in 

CM lO 0 



ID * 

0 

IO CM CM 

o 


►t ft 


• 

m • 

o 


cr * 




in 0 


b- ft 




O *0 


in * 

ii 

11 

II 11 

• • 


►ft « 




M3 

% 

o * 




<o 


« 

o 

LL 

I- U. 

O' 


UJ « 

o 

UJiO X 

4- 


x « 

»- 

x'k- o 

o 


3 ft 


c 

c z 

*4 


b- ft 



o 



< ft 



o 



cr # 


x: 



UJ * 


y- 




a. ft 


X o z 



y * 

UJ 

3 Z 

o * 



UJ * 

cr 

-j uj -* m 

in 


1 - * 

3 IX 3H Q 

—« 


« 


u -J 



z ft 

< H JUJO 

Lt 


-4 ft 

X 

< <,0 LL 

x in 


lx ft 

UJ 

u — 

UJ •-< 

■4 H 



CL 

X xx z 

»- 



X 

< < 

UJ 



UJI»- 

b- Z 

o o 



HZ HZ 

UJ 3 




< -• 

< Z 

in l- 



13*“* Z 

*-• o 

CL 



o,o 3,0 cr 

< 1— 



Zi< 

< u. 

_l -1 



urx cr oc 

UJ < 



cr 

UJ 

• 




•3 

• X 

• 0 




LJL- 

t- z 




UJ 

UJ 

o o 




X 

X 

r- o 


N 






Lmi 


Vo 


i *. 

15 

lb. 

U l 
;0 

u 
a : 
3 
b- 
< 
x 
!UJ 
X 


ILJ 

l> 

t • 

\b- 

l< 

[3 

13 

cr 


z 

o 


io 

ILJ 

icr 


* 

o 

!_i 


I- 

iMI 


l< 

ih- 

jl/) 

io 


IiONOKlI 

■ CO CT IT) 0 I 

■ r^- -ft cm cm| 
i i<r*cr o' a 


i0 0*0' dr 
® ^{< 


I iro 0-4 r- i 
I iCO 0 0 h* < 

:r- -ft: cm -* • 
I O' O'iO' O' ( 


|0 CM 0 < 

ico r-;-ft n r 


I Hz 

I Li CO 
| > uj< ►■« 


. - - !lj 

!-J 03 © 


iLi < 

a: x. 


J •! 
< H 
*-• in 
x *-• 

|< a 


to nl 3 n 
io cm m r- 








X 

i« 














I 

X 














X 

X 














3 

3 














»- 

h- 














0 

0 














0 

JO 














O 

O 














■f 

♦ 














r- 

—4 














0 

0 














0 

r- 














0 

r- 














r» 

-* 














CM 

-4 














• 

• 









a 















UJ 















3 





II 

II 









Z 




















‘H 

X 









t- 





X 

o 
















c 

_ 

c 

O' a 

rO CM 

0 3 ! 

c c 







>- 




o 

IO 0 

IO *ft 

0 3- 

O' 0 







J 2 

Z r-j 

X r- 

o 

O' CM 

0 O JV -ft. 

r- a 







H- O 

o u 

o o 

o 

3 ro 

r- -ft : cm *o 

0 0 


< 





z 

M4 

< 

1- 

o 

H IO 

* 0 

r» 0 



' — ' 







►- x 


• 

• • 

• • 

• • 

• • 


C3 





r\.j 

J L 













X UJ 

UJ 













3 

3 

3 








UJ 





u;u 

0 








X 


X 

ft 




X 








3 


uii 

ft 


»- 










CD 


>- 

* 


uJ! • 

>- UJ 










< 

» 


-J 

1- 

X 0 

H- 







LL. 



ft 


Z 

o 

3 3 

X 

O' '-ft o* 

0 r- 

0 0 

0 0 





ft 



b~ 

< b * 

X 

& 

0 H 

j’ in N a) 

O' CM 




z 

ft 




r 

r— x 

O' 

0 .0 

-4 CM 

-M 

© 0 





* 




w r— 

0 X 

r- 

0 H 

•ft 0 r- o' 

0 0 





ft 




Z -* 


O' 

o 13 

O' 0 ; IO O'ilO IO 




i- 

ft 




UJ z 

13 X 

O' 

O' o 

O' CM 

0 0 

O' 0 





* 




3 

3 3 

• 

• • 

• • 

• • 

• • 





ft 





>- 

CM 0 0 

CM 0 

O' 0 

-ft CM 




►- 

ft 




X 

0 


CM © 

O' 0 

0 0 

r- -ft 




o 

ft 




UJ 



•ft H 



1 




X 

* 




X 










X 

» 















ft 














o 

* 














*-« 

* 

■ r- 













o 

* 

b_ 













X ft 

4 



• 











* 

e 

L> 




0 

CM 

IO o 

,0 0 

r- 0 


in 

UJ:* 

0 

UJ 


ox 

X c r* 

ro ro 

p- r- 

CM —ft 



1- 

* 

X 

0 


UJ UJ 

0 2. 

C-.-* <M 




cm Cj cl — i!ro 

UJ: * 

U. 

X 


I- >- 


aio* O' 

O' O' 

O' O' 

0 O' 



ft 

cc, 

I- 


O «r 

UJ 

• 

• • 

• * 

• • 

• • 



• 


* 

X -J 

lL X 

cr 

a 3 










b. 



X 

X 

X 

X 









o 

« 

o o 

X 

X 











1 * 

r- o 

r* 3 

3 


3 









UJ: ft 

0 O 

0 »- 

1- 








ec rf'f- 0'ir 

cr 

ft 

• ft 

• tr 

0 


b~ 






IO 0 CM o 


3 

ft 

O' o 











ro 


* 

0 ^ 

010 


< 

0'CM 

0 © 

0 0 

P- 0 


O' 

< 

ft 

0 0 

O © 

3 


C^ft 

0 0 

r- p- 

CM -ft 


■ • 

• 

cr ! ■*» 

0 


3 

X 2 

r- 

-< CM 

—ft —ft 

0 © 

O' 0 




j * 



0 

< 

H- 

O' 

0* O' 

O' O' 

O' O' 

0 O' 




X 

* 


0 

O 

3C 

UJ 

» 

• • 

• • 

• • 

• • 




X 

« 


0 












Ul 

ft 


0 

O' 


X 






ro in!;* r- 

in 

K 

ft 


N 

0 








CO CM O h- 



ft 


CM 

H 


X 






cm cm cm o:ro 

UJ 

ft 



• 








O' O' 

O' O' 

O' 

Ulft 

ll II 

II 



13 






• • 

• • 

• 

< 

ft 














uT 

ft 


II 

II 

O 

b- 

o 

^ © 

0 -ft 

O' 0 

O CM 




X 

ft 




»— « 

X 

o 

•ft 0 

O 0 

CM -ft 

O' 0 




3 

ft 


o 


3 

r- 

o 

0 r* 

0 0 

3- cJ- 

CM CM 




0 

ft 


X 

UJ 

3 


o 

O' O' 

O' 0 

0 * O' 

O' O' 

CM CD >- CM 

0 

Ift 




X 


• 

• • 

• • 

• • 

• 

CD 0 Cf 4-4 


Q 

ft 

o o 

o 




•ft 






-* O 

IO 

Z 

ft 

o o 

o o 

>- 









O' 

< 

« 

>- X 

* 1- 

3 










• 


ft 


II 

r- 











0 

ft 


b- 

z 











UJ 

ft 

UJ 


UJ 


>- 










ft 

X 

b- 

X 


i- 









K ; ft 

3 

«■ <1 

X 



c 

cr 

r- cr 

0 ^ 

h* CM 

o 0 io r- 

0 

X 

ft 

r- ui 

>- 

3 


o 

o ( cm r- 

CM 0 

o- -ft 

r> 

io r-j r- cm -4 

UJ 

ft 

< X 

1- 

u 


o * 

©(O' 0 

co r- 

r- r* 

0 0 


3 9MiO 

X 

ft 

X z 

—4 



3 

oiO' O' 

O' O' 

O' O' 

0 * 0 * 


o 

ft lUJ 0 

u> 0! 


UJ 

• 

• • 

• • 

• • 

• • 

• • 



cr 

ft 

X 0 

O -111- 


> 

•ft 








CL 

ft 

X UJ 

3 Z'-* 












ft 

UJ X 

[UJ *H 

X 











o 

ft 

1- X 

> 

UJ 











*-• 

ft 














3 

ft 

UJ UJ 

UJX 



UJ 









J 

• 

u> u 

U UJ 



X 

o 

-4 -ft 

CM 0 

O- 0 

0 r- 

m it 

-4 J- 


b 

ft 

7T 2 

Z -A 


3 

©10 x 

o <\ 

a 0 

0 0 


0 O0 



UJ UJ 

UJ o 



0 

O O' O' 

O' O' 

O' 0 

0 ® 

0 f" 

r- O' 

*~4 


IX X 

X X 



0 X 

C3 

O' OV 

O' O' 

O' O' 

O' O' 

n 0iin in 

0 


IUI UJ 

UJ 



0 

o 

O' O' 

O' a* 

O' O' 

O' O' 

• • 

• • 




u u. 

X ►“ 



X 

• 

• • 

• • 

• • 

• • 






UJ Ul 

UJ Z 



X 

•ft 










X X 

X < 















3 















o 















o 











i b- 




c 










o n 








3 • 

© 

n o 

0 O 

n © 

n © 

o cM|in 

O X 







c h- 

o 

CM 0 

r- o 

CM 0 

r- 0 

in 0 r- « 

Q U 







►ft 0 

o 

-ft CM 

0 0 

0 r- 

0 0 

• • 

• • 

• — • 







X *-4 

• 

• • 

• • 

• • 

• • 










< a 





-4 


£ 



& 

IO 

n 


uj 

a 

in t 
X <c 

uj x 

’ ^ 
in a. 
x 

, ■ w 

leg H- 
3 

UJ . 

|Z *1 
M < 


a 

si 


3 

O 3 

3 * 

O' <7* 
•H O' 
ffi 


S-i 


un , 

i>H 

> • i 

, o. 

bs 

uZj uj! 

k >1 


|x tl 


SC ,o 

|o O' 

CM 

O' 

:* 

CO I 


in 

►h 

UJ 


. y 

O D| 
UJ 3 

in ♦ 

CL H 
< »- 
_J _J 
UJ <t 


2 « 

0 # 

i*-t « 
* 
CJ * 
UJ *, 
T » 
UJ # 
X 

"-s 

uj * 

1 *, 

*\ 

>- « 

r 

l< 


o c* 
in in 


cl co 
x xcc o d 


• IX 

> 


o id 

»- »-; . 

ffi m CD CQj 


ro -4 
o ol 
♦ 

O Ol 

o o 

O O ; 
-H CM*. 


* *, 
o o 

♦ ♦ 
•h m 

-4 O' 
0 O' 
3 & 


a q > c 
_i u7z a 
o azo a 
in 

a © ® in 


xu!l 
3 ol 


X 
UJ 

“■^*1 
x <!»-iuj 
< x|i~ o 

-J Xi< < 
O ZXI X 

in *-*;! O 

H- 

h- H *i/i 
z z z 
UJ ui>- > 
QOiQO 
►h *-*.o X 
o ox uj 
z zIuj z 

»-t *h!< UJ 







o 

u_ 

O 1 

r- o O' in 



X 





in i 

ffi) O 

o o 

# 


X 

CL CL 


UJ 

o 

r- i 

O' o';0' O' 

* 



I I 


X 

»- 

• i 

* • 

• • 

* 



X N. 


3 





* 


313 3 


H 

-J 




* 





<. 

<. 




# 

CL 

X 

0) ffi 


CL 

X 

o 1 








UJ 

X 

O 1 

-H ^ 

tX) tO 

* 

O 

i* 

in m 

!a 

O X 

C | 

r> in 

cm r- 


r- 

c 

o o 

!Z 

z 

O 1 

ffi) o 

•H O 

* 

o 

+ |4 4* 

IUJ 

UJ 

in i 

O' O' 

O' O' 

# 

• 

ffi) fO 



• i 

• • 

• • 

* 

O' 

=f 

K) CT 



o 




« 

in 

CM 

ro »n 


UJ 

z 




* 

ffi) 

3 

CM CM 


> 

< 




* 


• 

• • 


►H 


O 1 



# 





K 

in 

o | 

-h r-.o' cc 

* 





< 

i— i 

C 1 

r- ^io* »o 

* 

il 

II 

II II 


3 

o 

m i 

O -H 


* 





UJ 


CM 1 

O' O' i O' o* 

* 





a 


• f 

• • 

• • 

* 

o 

X 

Kt Um 







* 

o 

UJ 

O X 








1- 

ah a 







« 


0:0 Z 




C | 



* 



o 




o | 

*H -HI CM ffi) 

« 







O 1 

r*- o 

CM * 



Tj 





3 | 

3 -n 

rO CM 

« 


I- 





O | 

O' O' 

O' O' 

II 


X O’Z 




• | 

• • 

• • 

* 

LJ 

3 ZT 

c - 







* 

a 

_j uj;»- 4 in 







* 

z>iu. -i l- a 







* 

i- 


O -1 







« 

< 

H- _J UJ O 



z 




II 

X 

< < 3 U_ 


Ul < 

o 




* 

UJ 

UJ 

tj »-4 


> UJ 

»<H 


m cn 

ro 3 


0. 

X x,x z 


Xii- 

1 

•H CVJ 

O -H 


2 

< 

< 


*— 

u 


ro ro 1 r- ffi) 


UJV- 

H- X 


< * 

UJ 


r- 3 

3 m 


K'Z K 

4— 


-J 0 3 O 


• • 

• • 



< *H < X 


ul <;UJ 





UJ 

*+ z. *-• o 


X X 

X 





O O 3 0 X 








Zi< 

< u. 








ul 

XXX 








X 

ul 

— 








Ul 

• X 

• 3 


J • 



3 n 

3 n 


u. 

Ul 

♦- Z 


< 1- 



3 (M ffi) N 


UJ 

UJ 

o o 




O *H 

(M ro 


X 

X 

H- O 


X H 

NJ 


• • 

• * 






< o 






-4 3! it if 

i-« a . cc ^ 

%C O'' O O 
O' c0;O' O' 


\0 in ; ca 
* m- * O' cnj 

ooio C0-*-i 
O' O' 1 O' ffi : O' 


r- cm' co s-,-4 
cm O' *o in : cm 
— * O O O',*-* 
O' 0' i 0' ffi O' 


1 


=f- (M3 3i-4 
(\JM^ tOiCM 

OJ H H O H 

O' O'! O' O' [O' 


•| 


3 -H Q' ffi) CM 

in to o-i 3 
m in 'in n [in 


o 3:2 3 
o cm, m r- 
in ffi> ,r~ to 


3 *-4 





a 

X 










X 

X 










>* 

|\ 










3 

3 










1- 

►- 










®i 

CD 










in 

in 










o 

o 










■4 

4> 










r- 

—i 










3 

CM 










CO 

r- 










3 












o 










CM 

H 










• 

* 










II 

II 











r 










X 

o 












z 











c 

_ 

C* 

cv c 

a: 

CM CM 




>- 


►M 

< 

o;3 

<0 ffi 

3 3 




-1! Z 

Z l- 

X 

o 

CM 3 

3 r* 

3 ffi 




H»iO 

o o 

o o 

0:3 -« 

^ r- 

O' O' 




Z 


►H C 

H* 

o 

—4 fO 

cr 3 

3 r- 




u 

W 

►- X 


• 

» • 

• • 

• • 




g 

J 

J X 









X 

UJ 

UJ 









3 

3 

3 









(J 

UJ 

UJ 






cc 

» 



X 

X 






UJ 

* 










>* 

* 


X 

• 

>- UJ 






< 

* 


_l 

1- 

X ffi 

►- 






* 


ZlO 

-1 3 

X 

* 

o- 3 

3 ffi 

O' 3 


* 


M 

H- 

< r- 

X 

3' O' O 

-4 f'- 

3 O' 

z 

» 





)- X- ffi) O CNJ 

ffi 'M 

CM 4- 

o 

» 




r— v- 

O X 

3 3 0 

O 3 

-4 O' 

^H 

» 



iZ 

z — 

*H 

3 ffi 

3 3 

O* CM 

K 

* 




u z 

UJ V 

— H 

3 ^ 

3 r- 

0 r- 

VJ 

* 




3 

3 3 

• • 

• • 

• • 

UJ 

* 





»- 

r- cm *o 

3 ffi 

O' 3 

>- 

» 




X 

CD 

3,-4 O' 

ffi 

3 3 

o 

* 




UJ 



-4 



X 

» 




a 






X 

» 











* 










o 

* 










*“< 

» 

H- 









o 

* 

X 









X 

* 

4 



• 






CJ 

*• 

o 

CJ 


r- 



-H CM 

3 sr 

CM 3 

UJ 

* 

in 

UJ 


o a 

X 

f^ O CM 

cr cm 

r- * 

1- 

* 

\ 

in 


X X 

in x 

vO lO K1 

CM CM 

—4 -H 

uj: • 

X 

V 


»- >- 

H* 

O' . O' O' 

O' O' 

O' O'! 

2 

* 

cc. 



o < 

X 

• 

• • 

• • 

* *i 


* 

X -J 

X X 

X 

X _J 






u 

* 


X 

X 

X 

X 





o 

» 

Ip o 

*H \ 

\ 








* 

r*- o 

r- 3 

3 


3 





u 

* 

3 O 

CC >“ 

h 







cr 

4» 

• • 

• ffi 

<T 


X 





3 

* 

O' o 


in 







H- 

* 

in j* 

in 


< 

-4I-H CMJ 

3 a- 

CM 3 

<., * 

3 3 

o 

o 

-J 

M 

r*> 

O CM 

^ CM 

I s * * 

X' « 

to 

4- 

4 

-J 

X » 

3 10 3 

CM CM 

-4 *H 

UJI « 


r- 

O 

< 

H- 

O' . O' O' 

O' O' 

0* O' 

X 

* 


3 


X 

UJ 

• 

• • 

• • 

• • 


* 


a 

— H 







uu 

« 


3 

CM 


X 





»- * 



I s - 








* 


CM 



X 





UJ 1 » 


• 








o 

* 

II II 

II 



UJ 





<:# 










XI* 


II 

II 

o 

H 

o 

CM O 

O O 

3 3 

Xi * 





X 

0:3 CO 

r*- -h 

-4 P*- 

3:* 


O'M 

3 

h- 

0 0 

3 313 ^ 

in 

* 


XlUJ 

-J 


O O' O' 

O' O' O' O' 


« 




X 


• 

• • 

• • 

• • 

O'.* 

o o 

o 




H 




z 

* 

o o 

o a 

>■ 







<: * 

H X 

36 K 

-J 








« 


II 

t- 







im « 

• 

I- 

Z 







Ul * 

UJ 


Ul 


V 





LH 

« 

X 

►- 

X 


r- 





h* 

* 

3 - 

» <r 

a. 


►- 

C 0 CM 

C -4 

C 3 

X 

* 

K UJ 

V 

3 


u 

0‘3 O' 

3 -4, r* ^ 

UJ! * 

< X 

h- 

CJ 


o * 

O O' ffi 

ffi ffir~ r- 

X 


X 3 

•—4 h- 



J 

O O' O' 

O' O'! O' O' 

o 


uj in 

U X 



UJ 

• 

• • 

• • 

• • 

a 


x in 

O -J 

H- 


> 

•H 




n 

» 


_j rr 

V-4 









UJ X 

Ul 4-H 

X 







a 

N 

H x 

> 

UJ 


















3 


UJ UJ 

X X 



UJ 





3 


o o 




X 

0 

•M -H 

CM 3 

4- 3 

X 

* 

z z 




3 

0 

^ 3 

cr < v 0 to 


11U UJ 

Ul O 



n 

0 

O' O' 'O' O' O' ffi 



X X 

X X 



m Xi o O' o' 

O' O'] O' O' 


) UJ UJ 

UJ 



UJ 

O O' O' 

O' O' 

O' O' 



X X 

X 1- 



X 

• 

• • 

• • 

• • 



ul UJ 

uj Z 



X 

«H 






X X 

X < 











J 











o 











o 











CJ 














3 • 

3 3 3 

3 3 

3 2 







< r- 

O CM 3 

r- o cm 3 







*-t m 

O 

—4 CM 

3 3i3 r- 







X -H 



♦ • 

• • 







< Q 






.87 s .99866 .9712 .9386 ,9035 ,9035 71.472905 .9046(1 

.00U . 998 i >7 .9690 .9329 .9121 .9121 42.267320 l.OOQOn 

IEX1T1 



cr 






0 






IT 












•O 




1- h- 


r~ 




U. Ul 


K) 










O O 






UI Ul 






cr cr 






X X 

cr cr 

ui 




— ■ — * 

X X 

►H 




\ \ 

>*. Nk 





m J> 3 

3 3 

UJ 




I- H* 

y- >- 

a: 




CD CQ 

CD CD 

3 






i_n >- 




iO ^ 

* * 

a_ < 




O O 

0 0 

j a: 




+ + 

+ + 

V- UJ 




O 0 

C\| CM 

in cl 




a 0 

3* w 

z 




0 0 

-M 4- 

• UJ 




•H CM 

33 «M 

cc 1- 




• • 

• • 






| 

Ui * 






►- z 




II II 

11 II 

Z )"5 
















cc c 

> U 

M3 




-J Ul 

z cr 

O' 




0 .r 

0 0 





in *-« 

0 





0 0 

0 Ul 

% 










X 

X 

a 




3 

Ul 

UJ 




- _J 

jr * 

u: 




X Ll O Ui 

N 




3 

CL ►- 

h> 




-I Q 

< 

in uj 




u. ui'id cr] 

r~ 




cr 

z 1 

•O MJ 




cr <!MUJ 

r-i a 




< ari»— 0 

• •< 

• 



_l u.l< < 

vO O' 




0 ZibJ cr 

CM O* 




Ul Mil O 






H 

M3 




h- H* 

• Ul 

•H 




z z z 





UJ u 

3- 



2 * 


a a 0 t£> 



O * 


*m m. 0 a 



*-» * 


0 0 

cr ui 



H- * 


z z 

UJ z 



U» # 


►M »M 

< UJ 

m 


UJ * 




*-» 


"3 * 






UJ * 




Ui 


cr * 




z 


* 




»-• u 


1- * 




H- *- 


< * 






Ui « 




Ui >~ 


1 * 

* 


* 

> ^ 


* 






h- * 


K 


y- <j 


2 * 


U 


< 0 


< *, 




_i — 1 


»-» * 


cr 


Ui u. 


O * 


X 

cr cr 

cr > 


< # 


■ W 

X X 



cr « 


*s 

S X 



# 


3 

3 3 



0 « 


1- 

K H* 



z « 

cr 

CD 

CO CD 



< * 






« 

0 


in Ki 



z * 


o;o oi 

« » 


0 # 

•X) 


♦ + 



*-• # 

• 

M3 

0 0 

cr »- 


K ■» 

O' 


•M CO 

x u. 


3 * 

in 

CM 

W M3 



CD * 

M3 

tO 

CM CM 

0 


»-* « 



• • 

0 


cr « 




0 0 


K * 




r-4 O 


U) * 

II 

II 

II II 

• • 


-• # 




C 

, 

C* « 




a 


* 

0 

U. 

1- U. 

a 


UJ « 

0 

UJiO z 

a 


cr * 

>- 

crii- 0 

m 


z> * 


01c z 

»r 


H- # 



O 



< « 



O 



cr « 


X 




Ui « 


1- 




a. # 


x 0 z 1 



z * 

UJ 

3 Z 

O -1 



Ui « 

fj' 

-J UJ; »“* 1/1 

in 


»- « 

D LL J H Q! 

»-« 


» 

H- 


O -1 



z « 

<ih- -JiUJOl 

Ui 


-• « 

OCX <iD u-l 

z u- 


U- * 

UJI UJ *-« 

IU H 

►-» M 



Q.'IX,1£2 




Z! < 

< 

u 



Ui 1 )— 

I- X 

0 CJ 



hi2 HZ 

Ui X 




< 

< z 

n h- 



Ui 

►-* z;*m 0 

a. m 



OO D 

0 cc 

< H- 



Z|< 

< u. 

-J J 



ui cr cr 

cr 

ui < 



cr 

UJ 

• 




aJ 

• Q. 

• 3 




U. 

li. 

1- z 




UJ 

Ui 

0 0 




cr 

cr 

H- O 


» 

Lmi 


'u 

o 

tuj 

a: 

3 

«. 
cr 
u J 
a. 

s 


UJ 

cr 


< 

z 

,0: 

lo X 


ll/J 


I H z 
Ui < 0 

> Ul »-* 

U X,H- 

; u 
l< • Ul 
i-J C^O. 
ui <iui 1 
cr aiac 


j •! 

< *~| 

~ w , 

X M 

< a 



o o 
►« < 
f- *t 
J u 
UJ 


>- UJ, 

a. d 
-J 3f 
< »- 

15 


O Dll 
Ui UJ 

ss 


< 
LL »-« 

o o; 


X 

O X 
z — 

UJ x; 
-I 3 


CL 
UI X 


3*' 

UJ 

> 


UJ 

a: 

3 

lUl 

l/l SL| 
Ui 

a: 

0. 


-J • 
< ►— 
»-* l/>! 
X 


o c\ m 3 
o cm 
o o -Hj 
o in 
O'-* *oj 


^h 'o r- 
cm r- 

O XI o 

x) o 
po o* cm 

O' CO If) 

•’ • • 
K30* w 
O 


^-:r* ■*■ 
run vC. 
MJ CM cJ 
(MMT 


* 

ro in sO 

vD CM cv! 
o>:t* J* 


0 0 >0 
o in >o 

O CO N 

oio' a* 


•1 ■ 


C!»C IT, 
0 *0 CO 
O O' CO 
0|0' OH 


o!*h 

ciffl 
O O O' 
O O' J*| 
©O' O' 


n o 
CM m 
— < OJj 


* r* cd 
h- ^ m i 
^ ^!h ov 
if FVO* O'! 
& inivo H 


r* K) Ko oi 

H H’ CM f'*' 
vOM4- fMi 
co in O': 
10 0^0 
o r^jo r» 


m d 
f* of 

o d 

o o 
O' o 


O' CO 
o 10I 
m cm 

vO •<« 

ct- 

in co 


-.0 

cf no|o mj m 


M3 « 0 0|0> o 

«o m o* us. in ini 
h ►Joolp'C 

O' a* o* 0*1 - - 1 


M) lOliO O'O' O. 

<0 vn o' >o!cn m 

O' o 

CO O' 


O Oi 

O' O' 


-■ CM CM O' h* M3 

in co;<o cmiio r*- 
M3 lOJJt &ito cm! 
O' O' I O' O'O' O'] 


c « e *c;c M3 1 
* O* M3 CM: O* M3' 
X) r-|r- no m>| 
O' O'! O' O' O' o* 



<M lOl^ in,M>r*i 
=r <M- a to; .0 *[ 
O' O' ‘O* CO '33 CO 
O' O' | O' O', O' 0*1 

o' O' 1 0* O'; 0* o* 


X> ©[Ml 3|A3 

r- o 'cm in,N o 

*o injM> h- co o 




£ 



* 


X 






et 






O 






• 






r 4 




V- ft- 


* 




x ll. 


K> 










0 0 






in in 






X X 






X X 

X X 

IA 





X X 

•-» 




\ X 

\ \ 





3 3 

3 3 

Ui 




h- i- 

H 1 - 

x 




cn cd 

X CD 

' 3 
(A H- 




0 -« 

m in 

CL < 




0 0 

0 0 

LlI X 




♦ ♦ 

♦ ♦ 

I- UI 




0 0 

to to 

IA CL 




0 0 

r- co 

x 




0 0 

0 r- 

• Ui 




-« CM 

CM 

X K 




• • 

• • 

O 






u • 






ft- X 




II II 

11 it 

Z ft- 






*-« < 










a a 

> c 

£0 




3 UJ 

z X 

O 




0 x 

0 0 





in t~* 

UH 





0 0 

0 in 











X 

X 

O 




3 

UJ 

UI 




*■ 3 

^ **■ 

in 




X u. 

0 UJ 

V 




3 

Q.H 

H 




_i a 

< 

CO U. 




u. UI 

0 X 

H 




X 

2 

ft- co 

. 



X < 

*-« UJ 

vO <M 




< an 0 

• • 




3 IlK < 

CM CO 




O Z ;3 X 

tO ^ 




in 

xo 

in 





h- 

tO 




ft- i- 

• in 

CM 




2 22 





UU>>- 



z * 


OOOO 



0 * 


H M O QC 



*-4 * 


0 O X UI 



1 - * 


Z Z'uJ Z 



0 * 


M H 

< UI 

in 


UJ # 




*-• 


* 






UJ * 




UJ 


cr # 




z 


* # 




»-• in 


i- * 




1- — 


< # 






UJ * 




UJ v- 


X * 

* 



> »- 


* 




►— > 

ft- * 


ft- 


1 - 0 


z # 


IL 


< 0 


< * 




3 -1 


»-• * 


X; 

UJ UJ 


0 * 


xix a 

CL > 


< « 



X XI 



X « 


VN V 



« 


30 3 



c. # 


H'ft- K 



2 # 

X 

CCD CD 



< # 






# 

O 

:* in rO! 



2 « 

r- 

OOO 

«. » 


O « 

sO 

■ft 

4 * 4 * 



*-• # 

• 

>D CM 03 

CL ft- 


ft- # 

O' 

jfiw in 

X u. 


3 « 

in 

n:h 0 



cr # 

>0 

IO.CM CM 

0 


# 


• 

• • 

in 


a * 




CM CM 


ft- # 




^ ro 


*n # 

11 

ll 

II II 



*-« « 




X 


0 * 



» 

0 


* 

0 

tL 

1- u 

m 


UI # 

0 

u,o z 

fO 


X « 

►— 

X: 4 — O 



3 * 


CIC 2 

ro 


ft- * 



0 



< « 



0 



X 11 


X 




UJ « 


ft- 




a. * 


X IS 

Z 1 



z * 

UJ 3 2:0 - 



UJ * 

X 

-J uj’*-« in 

in 


ft- « 

3 

u 3 ft- 0 

H 


« 

ft- 


0 -J 



2 * 

< 

»- _J l*J 0 

UJ 


<-* « 

X 

< <ou. 

z >n 


Ll * 

UJ' tu *-« 

UJ *-« 

-« M 



X'l XlX Z 

H 



2 

<z 

< 

Ul 



UJ , ft- 

>- z 

c 0 



ft- 

2 ft- 

z 

UJ 3 




< *-*!< Z 

m r- 



UJ 

H 2 

-• 0 

CL 



0 

0 3:0 X 

< ft- 



2 

< 

< u. 

-J -1 



UJ 

X XX 

UJ < 



X 

UJ 

• 




3 

• X 

• 3 




U 

IL 

H 2 




UJ.UJ 

O O 




X 

X 

H- O 


<M 

Q ( • 


1 

' 











X 

X 


















X 

X 


















\ 

N 


















3 

3 


















ft- 

ft- 










. 








X 

03 


















in 

in 


















0 

0 


















♦ 

+ 


















r- 

0 


















vO 

CM 


















CD 

O 


















vO 

O 


















r- 

CM 


















CM 

-4 


















• 

• 












c 



















U. 






11 

II 












2 



















z 






H 

X 


















' X 

L 3 




















Z 











z 








c 

- 

c 

— 0 

cv 0 







G 






>- 


»-* 

< 

O' sO 0 

O' ft- 







L- 






3 

z 

2 ft- 

U_ ft- 

O O' ft- 

O' 4 ft 













t- 

O 

O <J 

0 0 

o? a 

to xO 













Z 

*-» 

►H < 

ft- 

0 

«H IO 

in 







<1 






0 - 

— X 


• 

• • 

• • 







'w 







4 

4 L 











C 3 






a uj 

3 

















3 


> 

















U UJ 

3 














X 

* 



X 

X 











X 



UI 

# 


»-l 















>• 

# 


UJI • 

>- UJ 











c 



< 

# 


3 ift- 

X 23 

ft- 













3 

* 


ZIO 

3 5 

X 

m cM * 

(M 







uu 




# 


*-• 

ft- 

< ft- 

X 

•H H ft 

CM O' 










Z 

* 



IX 

ft- Xi :m ft- C\J 

ft- M) 










O 

* 




ft- ft- 

O X 


xO M 3 










t -4 





2 -4 

z - 

O' O' 0 

rO O 










h 

* 




3 2 

3 V. 

ft- O ^ 

X M> 










O 

# 




3 

3 3 

• 

• • 











UI 

# 





ft- 

-O ^ ft- 

O' CM 










►— 

* 




X 

CQ 

<0 

CM 0 

O* O' 










O 

* 




UJ 














X 

» 




a 














X 

* 



















L* 


















0 

• 


















M 

# 

ft* 









z 


4 






0 

* 

u 









0 








X 

* 

« 



• 






H 

O 







0 

* 

0 

O 




ftO 

in cm 

*H CO 


►* 

O 


to >c 

in ir 

'OhWC' 

H 

UJ 

• 

in 

UJ 


0 X 

X 

CO CM Cd 

^ O', 


O 

O 


CO CM O' d- 

cm O ' 1 >o in 

in 

y~ 

* 

V 

in 


u u 

in s 

in cm ci 

-H O 


UI 

O 


in O' ‘O' O' 

O' cc cc o' 

LJ 

* 

u 

N. 


ft- >- 

ft- 

cr 

O' O' 

O' O' 


X 

O 


O' cc <C cc 

CC X 

co coico 

z 

ft 

X 

1 - 


0 < 

3 

»■ 

• • 

• • 



• 


• w 

• « 

• • 


• 


* 

X 3 

U X 

X 

X 3 





2 

0 

*-1 







UI * 


i 

X 

X 

X 




H 









O' * 

O O 

H \ 

\ 






u 










* 

r» 0 

ft- 3 

3 


3 





c 

c 








¥ 

c. 

X ft- 

ft- 






u 

3 

0 

inccioo' 

0 O'; 3 J- ^ 


X 

ft 

• • 

• ffilCD 


ft- 




0 

U. 

0 

1 

CO ^ 

in 

cf o r» r- : m 

3 

ft 

O' 0 








UI 


in 

1 in O' 0 o' 

O' O' CC f'-’O' 

H 1 * 

in cr 

inin 


< 

ro 

in co 

•H to 

O 

r~ 


O' <0 O' CC 

CO cc co CO CO 

< 

ft 

M 3 M 3 

oio 

3 

»-• 

X 

CM CM 

* 0 

X 

H 

• 

1 

• • 

• 

• • 

• • 

• 

am 

cc 

-f 

-ft 

3 

X X 

in CM CM 

•H 0 

3 



1 






al 

ft 


r- 

•HI 

< 

ft- 

O'lO' O' 

> 0 

H 

3 


1 






a. 

ft 


sO 

eft 

X 

3 

• 

• • 

• • 

< 

< 


1 






z 

ft 


a 

cc 






X 

Z 

0 

1 






UI 

ft 


M 3 

vO 


X 




UI 

X 

0 

1 

ro -4 

CM 

O' -4 

in m -M 

K 

ft 


ft- 

in 






X 

0 X 

0 

1 

CO 

vO O 

r- ^ 

0 o ! m 


ft 


CM 

H 


z 




Z : 2 

c 

1 

in a 

0 c 

O' O'* O' cc 

O' 

UJ 

ft 


• 

• 






UJ 


n 

1 O' O' 

O' O' 

•O 30 

3 co 

< 0 

u 

ft 

II II 

II 



3 




l- 

UI 

• 



• * 

• • 


• 

< 

ft 










0 








u 

ft 


II 

II 

Q 

ft- 

0 

ft- in 

ft- 0 

UJ 

2 








a 

ft 




*— < 

U 

0:^ in 

ro in 

> 

< 


1 






3 

ft 


O 

*-« 

3 

ft- 

0 X ft- 

<0 in 

H 


0 







m 

# 


XIUI 

3 


0 , 0 ' O' 

0 0 

>- 

in 

0 


to co sO 

lO CO 

O' O' 



ft 




U 


• 

• • 

• • 

«I 

►-« 

0 


cr 0 k> \0 

fO CD 3 t 3 T 

in 

01* 

0 0 

0 




•H 



-1 

0 

_n 



—l C 

0 O' O' to 

O' 

z 

ft 

0 0 

0 0 

>- 






u 


CM 


O' O' O' O' 

O' cO;CO CO 


< 

ft 

ft- X 

3 K ft- 

3 






0 


• 


• • 

• • 

• • 

• • 

• 


ft 


II 

K 















in 

ft 

» 

ft- 

Z 









1 






UJ 

ft 

UJ 


3 


>- 













P -4 

ft 

X 

ft- 

X 


ft- 






c 







t- 

ft 


* -IIX 



c 

^ c 

CM ft-. 



0 


m in cm ^ 

co in cm 

H 

X 

ft 

P UJ 

>- 13 


O 

0 

to x 

to X 



0 


-o CM 

cm cr 

O' n O' O' n 

UI 

« 

< X 


<3 


O * 

O O' co 

x r- 



0 


n cm cm — 

0 oio' «!o» 

X 

ft 

X 3 

ft- 



3 

0 O' O' 

0 0 



0 


O' O' CT' O' 

O' O' CO cC cO 

0 

ft 

ui in 

O UJ 



3 

• 

• • 

• • 











* 


vl. ~~t 

ft- 


-> 













X 

ft 

>■ Ui 

3 2 

H 
















ft 

UJ X 

UJ »-• 

X 















a 

ft 

ft- X 

J> 

3 















*-« 

ft 










z 








3 

« >UJ UI 

UI X 



3 




UJ 0 








3 

* 

0 3 

a uj 



X 

0 

—4 — « 

CM O 

> UJ *-« 


1 

rO ^ CC CN 

^ r- 

IO CM 


u 

ft 


2 k 



3 

p X X 

^ CM 




CM O' CO J 3 

^ O' 

10 in 

n 



3 UJ 

3 0 


in 

O ' O' O' 

O' 0 

H- 

3 



O 3 

3 - K 1 

CM 0 

0 0 

»-4 



X X 

X X 



in x 

O O' O' 

O O 

<1 • 

3 



r*» o.m in 

mm n in m 


!UJ UJ 

>3 



3 

O! 0 ' O' 

O' 0 

3 0 3 0 



• • 

• • 

• • 

• • 

» 


.u u 

U ft- 



a 

• 

• • 

* • 

UJ < UJ 


1 








UI UJ 

3 2 



X 

-H 



X X 

X 


1 








oc <r 

X < 










1 









3 



















pi 










1 









0 










1 









a 







3 • 




3 n 

0 n 

o n 

3 ,n 

3 -• 







3 • 

0 

n o 

n 3 

< »- 




O CM 

in 

O CM 

in i*- 

O X 







< ►- 

0 cm m 

r- 0 

•h in 




O -« 

CM to 

m >0 

r- co 

0 UI 







~ in 


-4 CM 


X M 

•M 



• • 

m • 

• • 

• • 

• *>» 









• • 

• • 

< O 















< 3 





5 


■ M> 
t 30 
I eft 


« CM 
'O' 
'CO 

> c 


) in H 
' m O' 
jIo o»' 

MO* <0; 


1 <0 

'O' 

>;cm 

MG' 


t ) 


) ) o o 


3 


) « 


.000 .99647 .964? .9229 .6951 .6951 55.766302 l.OOQOn 

«EX1T) 



I ( t 


* 






K> 






O' 






• 






f" 




►- ►- 


IQ 




U. Ul 


tO 










© o 






in in 






X X 






X X 

X X 

U) 




«-* w 

X X 

H 




X x 






X X 

X X 

UJ 




V- h- 

►- y- 

x 




CD (X 

(D CD 

x 






in *- 




to ** 

in in 

a < 




o o 

o o 

J X 




♦ + ! + ♦ 

V- UJ 




o o 

r- in 

in a. 




O 0 (0 IQ 

x 




o o 

O -4 

• ui 




-4 CM 

to CM 

X K 




• • 

• • 

O 






UI • 






h- z 




II II 

II II 

z »- 






M < 










X c 

> IS 

o 




J UJ| 

Z X 

•H 




O X 

o o 

H 




in »-» 

o y- 





o c ? 

<3 in 











X 

X 

u 




X 

UJ 

UI 




- -1 

M - 

in 




X u. 

O UJ 

X 




X 

cl t— 

h 




-1 o 

< 

O' U. 




U. uilO X 

in 




X 

z 

O' o 




X < 

*-4 ui 

-4 O 




< a t- o 

• • 




-1 u. 

< < 

to 3 




o z 

Ul cc 

*Q fO 




in ►-* 

I o 

CM 





►- 

3 




H- ►- 

• in 

CM 




Z Z!Z 





uj uji >- H 



z * 


X X 

a <d\ 



o * 


M M 

o X 



►4 * 


O O X UI 



H- * 


Z Z ul Z 



u * 


»-* K4 

< UJ 

in 


UJ * 




*-» 


"5 * 






UJ ft 




UJ 


cr * 




z 


* 




*4 in 


H ft 




H- *-• 


< * 






X ft 




UJ >- 

• 

X * 




> 1- 


* 




►4 *-* 


H * 


K 


t- u 


Z * 


U 


< o 


< ft 




_l -1 


« 


X 


UJ UJ 


o * 


X 

X X 

a > 


< ft 



X X 



X ft 



\ % 



* 


x:x x 



c * 


K’K K 



z * 

x 

COS CD CD 



< ft 






ft 

o 

if 

in to 



z ft 

r- 

o 

o o 



o * 

® 

♦ 

* + 


l~ ft 


® 

(O H 

x k 

>- * 

O' 

3 

CM -4 

x u. 


z> ft 

in 

CM 

^4 O' 



X ft 

® 

to I CM •-» 

o 


*4 * 


• 

• • 

r* 


x * 




<\i o 


»■ « 




~4 o 


in » 

M 

II 

II II 

• • 


*4 « 




o 


c ft 




CO 


« 

o 

X 

►- X 

O' 


UJ * 

o 

UJ 

o z 

o 


x « 

►- 

X 

H- Q 

*-> 


X * 


c 

►C Z 

to 


»- ft 



O 



< * 


• 

o 



a * 


X 




UJ ft 


h 




a. # 


x o 

z 



Z ft 

UJ 

x z 

o * 



UI ft 

cr 

J UJ 

»-» in 

in 


►- « 

X 

X -1 

K O 

►4 

ft 

K 


U -J| 



Z ft 

<:»- JUO 

UJ 


-4 ft 

cr !< < 

X X 

z in 


U. ft 

UJ: UI *- 

x *-• 

-4 M 



CL 

X x 

X z 

►" 



z 

< 

< 

UJ 



UI’V- 

*- z 

o a 



HlZK 

z 

UJ 3 




< *-« 

< z 

in h* 



UI 

►4 2 

M O 

0. ►-» 



UlO 3 0 X 

< H- 



z ;< 

< UJ 

-1 -I 



uiicr xix 

UJ < 



x 

Ul 

• 




j 

• X 

• X 




in 

X 

y- z 




UJ 

UJ 

o o 




X 

X 

y- o 


? L/t7TT 


K- 


Ir 


!0 

I 

X 

jx 

p 


15 

a 

lui 




*< 

:_j 

Ui 

x 


o 


° 
I m 


x 


o 

o 

o 

o 

in 


o 

o 

c 

m 

CM 


C 

o 

o 

o 

o 


UJ <1 


_J C3 
UI <1 

x xi 


z 

o 


o 

UJ 

X o 
UI ! 

0 c 


-I 

i< ►— 
•-« in 
x *-• 
i< O 











X 

X 
















Zll 
















\\ 
















x^x 
















1- 

















CD' CO 
















inta 
















oo 
















3 

















r- 

c 
















® -4 
















<c 

CM 
















o 

3 
















n- 

CM 
















CM 

—4 
















• 

* 
















II 

<1 
















*-4 

r 











c 





I 

Q 

















z 

















c 

_ 

c. 

to O 

U. m cm k 

3 C 









>- 


►4. 

< 

O DO h* 

3 «-4 

cm in © o 









_J z 

L_Z *— 


3, O' £ 

O' X 

3 n »o o 









►-o 

O <J 

O O o 

3 O 

*0 3® O' O o 









z 


I-* < 

y- 

o 

-4 IO 

3 in.® f- 

O' o 









lli- 

*- X 


• 

• • 

• • 

• • 

• • 




c 





X 

J 

J X 















X-UJ 

UJ 
















■> 

3 















uw 

UJ 












X 

ft 


lx 

X 










c: 


Ul 

ft 

















ft 


UJ 


>- UJ 















-J 

H* 

x ir 

*- 











_l 

ft 


z 

O 

J X 

X 

tO 

3 o 

o ® 

rO tO 

o o 









M 

1- 

< H 

I 

in 

tO 3 

h- 3 3 ® 

3 *4 






z 

ft 




X 

*— x 

c 

CM O' 

® ® 3 *0 

r- —4 






o 

ft 




i- y~ 

3 X 

to 

3 r*> 

® n r- n 

•H O 











-■r m 

z — 

r- r- r* 

r- n 1 -* o 

CM —4 






1“ 

ft 




u z 

UJ \ ^4 O' -4 

o r- r- co 

h* in 







ft 




X 

-IX 


• • 

• • 

• • 

• • 













in3N 

O' CM 

3 o >4 r- 






y- 

ft 




X 

CD 


CM o 

O' O' 

® ® 

® in 






o 

ft 




X 



w-t H 









X 

ft 




a 












X 

ft 

















ft 
















Q 

ft 
















•-4 

ft 

y- 















o 

* 

X 
















ft 

4 



• 












o 

ft 

© 

u 




vC vO CM 

o r» 

3 O' 

»4 H 


O' r- 

O' O' 1 3- CM ^4 

UJ 

* 

l/l 

UJ 


o X 

X 

f* 

CM CM 

3 O' 

K> ® O' 3 




3 


ft 

V 

in 


UJ u 

in z 

in CM CM 


© O' ® o* 

in a 

O O' 

O'O'CMT 

Ui 

ft 

X 

\ 


K > 

H 

O' O' O' 

o O' 

O' CO 

03 C 

O'C-O'CC 

CC tz.to 

Zlft 

cc 

I- 


o c 

X 

• 





• • 

• 4 

• 1 

• « 

• 


ft 

X -J 

X XX 

X -J 












IL 

ft 


II 

X 

X 











O 

ft 

o o 

w N. N. 














ft 

r- o 

X X 


X 











UJ 

ft 

M3 O 

cc »— 1»— 








vO KJ 

xO CM 

CM O 

in to 


X 

ft 

• • 

• to:® 


i- 






r- r-iro co 

sO IO O' O' 

3 

X 

ft 

O' © 










in O' 

O O' 

O' O' ee r- o 

y- 

ft 

in 3 

in m 


< 

vC vO CM 

o r- 

3 O' 

•4 »4 

0 X 

0- U, 





4* 

v£ £ 

c 

© 

-J 

1-4 

r- 

CM CM 3 O' 

*o cc 

O' 3 





• 

X 

ft 

CC 

3 

♦ 

-1 

X 7 

in CM CM 

-4 o 

O O' 

® O' 






aj 

ft 



O' 


►— 

0* 

O' O' 

O' O' 

O' ® 

® ® 






CL 

ft 


sD 

3 

-* 

X 

• 

• • 

• • 

• • 

• • 






z 

ft 


a 

I® 













UJ 

ft 





X 






vO CM 1 »-4 CM 

co o *o to 


h- 

ft 



in 








h“ IO CO CM 

O' sO CM CMi 3 


ft 


CM 

^4 


z 






in o . o o 

o' o'’ O' tfi: O' 

Ul 

ft 


• 

• 








O' O' 

O' O' 

co X> ® 

C 

o 

ft 

II II 

II 



X 








• • 



< 

ft 
















IL, 

ft 


II 

II 

o 

H 

oMn 

® ® 

in cm 

r- r- 






X 

ft 




►- 

X 

o 3 in 

to in in O' 

O' CM 






X 

ft 


Oil-* 

X 

t- 

Oi® r- 

® in 3 in 

CM CM 






in 

ft 


ZIUJ 

-1 


OIO' O' 

O' O' 

O' O' 

O' O' 

vO CM 

O' **4 

O' o 

-* CM 

-4 


ft 




X 


■ 

• • 

• • 

• • 

• • 

r- cm 

^ cc 

t o »C vO 

3 

o 

ft 

© o 

o 




•H 






*■4 o 

3 O O' O' O' 

"Z 

ft 

o o 

o o 

>- 








O' O' 

O' O' 

O' O'.® ® 

® 

< 

ft 

1- X 

* M-l 











• • 



ft 


It 

:h- 













in 

ft 


K 

z 













UJ 

ft 

Ul 


UJ 


V 











*-4 

II 

X 

l~!X 


H- 











HlK 


- <ia 


*— > 

o 

3 C 

cm r> 

® cc 

O *4 

vO vO 

CM O 

N 3 0' ^ 

•-< 

X 

ft 

y- uj 


X 


U 

o 

*o ® 

*0 C0‘3 O 


h- cm 

CM ^ 

O' (O CO O' 

3 

UJ 

ft 

< X 

»- 

U 


o * 

O O' CO 

® r- r- r- 

vO 

n cm 

CM ■H 

3 0 0' C O' 

X 

ft 

X X 



-J 

O ; O' O' 

O' O' 

O' O' 

O O' 

O' O' 

O' O' 

O' O' ® CO ® 

o 

ft iuj in 

o UJ 



X 

• 

• • 

• • 

• • 

• • 

• • 


• « 

• • 

• 

X 

« 

x in 

O -III- 


> 

*4 










X 

ft 

Z Ul 

-1 Zjftft 














ft 

Ul X 

Ul M X 













o 

ft 

I- X 

> 

UJ 














« 
















X 

ft 

UJbj 

U X 



X 











-1 

« 

L> U 

a uj 



X 

O 

-4 -H 

CM *0 

* n 

® r» 

-4 in 

It «£ 

c in cm c 


u 

« 


Z J 



X 

©i ® m: 

3 <\ 

© cc 

'C 3 


(O O 

® C 

CM 



UJ UJ 

UJ o 



in 

O |0» O' 

O' O' 

O' ® 

® ® 


rg ho c 




X X 

X X 



in x 

O O' O' 

O' O' 

O' O' 

O' O' 

r- n in in 

in in n in 

m 



U UJ 

Ul 



UJ 

o 

O' O' 

O' O' 

O' O' 

O' O' 

• • 

• • 

• • 

• • 

• 



X X 

X K 



X 

• 

• •: 

• • 

• • 

• • 








UJ UJ 

UJ z 



X 

-4 












X X 

X < 

















-1 

















o 













- 




o 













1- 




u 










3 n 

o n 

o n 

3 -• 







J • 

3 

•n 3 

•n 3 

® 3 

n o 

O (M 

uo r- 

O CM 

n r- 

O X 







< K 

3 

cm in 

r- o cm a. 

r- o 

O »4 

CM to 

in sc 

h> ® 

O UJ 







*-4 m o h(\ 

to ® 

® r-i® o 

■ • 

• • 

• • 

• • 

• M* 







X «-* 

• 

• • 

• • 

• • 

• • 





•H 







< X 





*4 


-4 

X 

UJ 


cn 


• i > > j ) 


■ i 




line 



52 


b. Reentry 

(i) Mission Description 

The reentry altitude and velocity profiles used for the system 

output shown in Figures 13 and 14 are given in Figure 12. 

Input parameters for reentry are identical to those given in the mission 

description for ascent except for IC0NV which is set equal to 2 to 

designate reentry. The system parameters for reentry are printed prior 
to the current system output shown in Figure 13a. 

(ii) System Response 

The response of the radiator system during reentry is summarized 
in Figures 13 and 14. Figure 13 (a) provides selected program output 

i > , t 

showing current system conditions at various times during the reentry 
phase of the mission. Figure 13 (b) is a plot of the total heat 
rejection from the radiator system and the temperature of the coolant 
fluid at the exit plane of the flow channel. Figure 14 is a graph 
of the aerodynamic heating of the system as a function of time during 
the reentry phase of the mission. 
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c. Orbital Mission 

(i) Mission Description 

Transient system performance is predicted in Figure 16 for a typical 

space base orbital mission. The mission is described by a 270 nautical mile 
* 

circular orbit inclined at 55°. Incident flux values for solar irradiation, 
earth albedo and earth irradiation were supplied by a previously executed MRI 
program with the flux data being transfered from magnetic tape. 

Incident radiant flux values for a typical fin element are plotted 
as a function of time into orbit in Figure 15. Other input parameters which 
specify system geometry, material properties and fluid flow conditions are 
identical to those used for the ascent and reentry phases that were previously 
reported in sections 5a and 5b . 

(ii) System Response 

The response of the radiator system during the first twenty minutes 
of the orbit is summarized in Figures 16(a) and 16(b). Input parameters 
and a selected number of output listings of current radiator conditions are 
shown in Figure 16(a). Figure 16(b) graphically summarizes the outlet coolant 
fluid temperature and the total net heat rejection from the radiator for the 
same period during orbit. 
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Fig. 16(b) System Response For Orbital Mission-Heat Rejection 



d. Steady State Environment 


(i) Description of Environment 

Program output has been produced for steady state system performance 
to a time invariant environmental conditions. The entire radiator system is 
initially assumed to be at a uniform temperature of 200°F (TI = 659. 67R) . 

The integration process is started and it proceeds until steady state operation 
of the system is achieved. The environmental conditions during the integra- 
tion process remain constant. For the output of Figure 17, incident radiant 

flux values were read in from cards according to the format described in 

2 

section A.3(d). These flux values were 100 Btu/hr. ft. for solar irradia- 

2 

tion, 2.0 Btu/hr. ft. for earth albedo. There was assumed to be no infrared 
earth irradiation. There was no aerodynamic heating of the radiator panel. 

All other input parameters are the same as those used for program output 
during ascent (see section A.5(a)). 

(ii) System Response 

The response of the radiator system during the period leading to 
steady state operation is shown in Figure 17 . Figure 17 (a) shows the input 
parameters and a selected number of output listings of the system performance. 
Figure 17 (b) graphically summarizes the resulting outlet coolant fluid tem- 
perature and the total heat rejection of the radiator system as the system 
approaches steady state operation. 
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17(b) System Response for Steady State Environment— System Heat 
Rejection And Outlet Fluid Temperature 


B. Program Description 


1. Structure 

The program structure is dictated by the calling sequence of the 
Runge-Kutta-Simpson integration procedure (see Section III-2 of the System 
Analysis Manual). The Source Deck consists of two parts. The first part 
is referred to as the Permanent Part and the second part is called the 
Selective Part. The Input Data Set follows the Selective Part. 

a. The Permanent Part 

The Permanent Part requires no attention from the user unless he 
plans major program modifications such as consideration of geometrical 
configurations other than the one considered on the present analysis. The 
Permanent Part contains the first 42 program units. 

b. The Selective Part 

The Selective part contains the options for selection of system 
materials. It is assembled from case to case by the user in accordance with 
his choice of the coolant fluid material, flow channel material, meteoroid 
protection layer material, fin panel material and thermal control coating 
material. There are always 20 program units in the Selective Part. 

c. The Input Data Set 

The Input Data Set follows immediately after the Source Deck. 

As discussed in Section A-3 the Data Set consists of six groups that specify 
the non-symmetrical panel conditions, ascent and reentry profiles, the 
irradiation on the panel, inlet coolant fluid conditions, system specifications 
and finally the program options. These separate data units are arranged in 
an order of increasing possibility of input changes so that data units most 
likely to be changed are at the end of the Data Set. 



2. Fastrand Execution 


A run procedure was written for the purpose of entering the entire 
program, including property subprograms, on FASTRAND drums and collecting, 
or mapping, the program into convenient elements. The elements are chosen 
so that they form commonly used groups of subprograms which then may be 
selectively collected into a complete program. The advantages of this 
technique are the significant reduction in the compilation and collection 
time and the elimination of submitting the program in card form. 

The entire program is mapped into two basic elements: the perma- 

nent element and the property elements. The permanent element (called PERM), 
consists of all subroutines listed in the permanent part of the program; 
the property elements consist of the material properties of the surface 
coating, fin material, tube material, meteoroid protection material, and coolant 
fluid • 

Each property element has a unique name. The first two letters 
of the name are a prefix which indicates the function of the particular 
material. For example, a property element starting with the prefix CF 
indicates that those properties to be used for the coolant fluid. The 
remaining letters of the element name are a suffix which describes the 
particular material to be used. For example, the suffix HE indicates the 
material helium. The prefix-suffix combination CFHE, therefore, denotes 
the coolant fluid of helium. 

The sixteen names designating the permanent part of the program 
and the property elements are as follows: 

i) Permanent Elements 

1. Element PERM contains all subroutines in the permanent part 
of the program. (See Part B.3 a through g) 

ii) Property Elements 

A. Coolant Fluid Properties (Prefix CF) 

2. Element CFHE contains all helium properties 

3. Element CFSIL contains all silicon oil properties 

4. Element CFNAK contains all NaK properties 
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5. Element CFFC43 contains all properties of the 3-M fluoro- 
chemical liquid FC43 

6. Element CFFC75 contains all properties of the 3-M fluoro- 
chemical liquid FC75 

B. Fin Properties (Prefix FN) 

7. Element FNAL contains all aluminum properties 

8. Element FNCU contains all copper properties 

9. Element FNBR contains all beryllium properties 

C. Tube Properties (Prefix TB) 

10. Element TBAL contains all aluminum properties 

11. Element TBCU contains all copper properties 

12. Element TBBR contains all beryllium properties 

D. Meteoroid Protection Properties (Prefix MP) 

13. Element MPAL contains all aluminum properties 

14. Element MPCU contains all copper properties 

15. Element MPBR contains all beryllium properties 

E. Surface Coating Properties (Prefix SC) 

16. Element SCZ93 contains all zinc oxide coating (Z93) properties 

a. Control Statements 

The run deck consists of the following control statements: 

The catalogue card 
@ CAT ,P NASA. ,F2 

is used to create a file, called here NASA. This card is to be used 
only the first time the FASTRAND program is run. It must be omitted 
when subsequent runs are made because the file NASA should be created 
only once. 

The assign card 
@ ASG,AX NASA. 

is used to name an external file which has been previously catalogued 
and to cause its assignment to the run. The option A specifies that 
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the file is currently catalogued and the option X specifies that the run 
is to have exclusive use of the file until the run terminates so that no 
other run may interfere with the program operation. 

The use card 
@ USE N., NASA. 

allows the user to refer to the previously catalogued file called NASA 
by the single letter N. The abbreviated file name will save a great 
deal of card punching for the subsequent portions of the run. 

The Fortran compiler cards, one of which would be 
@ F0R,IS N.MAIN 

creates Fortran symbolic and relocatable elements both called MAIN 
which are to be entered into the file NASA. The I option causes the 
insertion of the element MAIN and the S option produces a listing of 
the element. There must be a F0R card for each subprogram that is to 
be entered into the file NASA. 

After all of the subprograms have been filed, they are mapped into separate 
elements. This is achieved by the MAP and IN card. 

The MAP cards, one of which would be 
@ MAP , ISR ,N.PERM 

specifies that the collector should combine a set of relocatable 
elements from the file NASA into a single relocatable element called 
PERM. There must be a MAP card for each program element that the 
user expects to utilize during program operation. A listing of 
these seventeen elements is given at the beginning of this section. 

Following each MAP card is a list of the separate subprograms that are to 
be placed in the element declared by the MAP card. The placing of each sub- 
program into the element is achieved by an IN card. 

The IN card, one of which would be 
IN N.MAIN 

places the element MAIN from the file NASA in the element name 
declared by the previous MAP statement. There must be an IN card 
for each program element that is to be entered into the element 
declared by the previous MAP statement. 
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As an example, if the user wishes to collect the subprograms called AA, BB 
and CC which are located in the file NASA and place them in an element 
called ELEM, the proper cards would be: 

@ MAP , ISR , NASA. ELEM 
IN NASA.AA 
IN NASA.BB 
IN NASA. CC 

It is necessary to place a define card at the end of each mapped element. 

A define card, one of which could be 
DEF THCFN , DTHCFN , CPFN 

is used to list those external definitions to be retained by the 
resulting relocatable element. The DEF statement causes the collector 
to construct a table defining the entry points to the element. 

There should be a DEF statement following each MAP card which 
declares the entry point to each element collected by the MAP 
card. 

b . Input Data 

The runstream described above enters all subprograms into the file 
called NASA and then subdivides the program into convenient groupings of 
elements. Once this run has been made, it becomes a relatively simple 
task to collect and run a typical system simulation. As an example, 
assume the user wishes output for a radiator system consisting of an 
aluminum fin, copper tube, beryllium meteoroid protection layer, silicon 
oil coolant fluid and Z93 surface coating. Assuming the user has pre- 
viously declared the file NASA and he has mapped the sixteen elements 
discussed above, then this combination of materials could be mapped into 
an absolute element called VER1 by the runstream consisting of 

(Run cards appropriate to the users system) 

@ ASG,AX NASA. 

@ USE N. ,NASA. 

@ MAP, IS ,N . VER1 
IN N.PERM 
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IN N.FNAL 
IN N.TBCU 
IN N.MPBR 
IN N.CFSIL 
IN N.SCZ93 
END 

@ XQT N.VER1 

This deck must be followed by the data set described in 
Section A-3. 

As a second example, if the user wishes to run the program for a,,., 
system consisting of an aluminum tube, fin and protection layer with helium 
coolant fluid and Z93 surface coating the deck would consist of 

(Run cards appropriate to the users system) 

@ ASG,AX NASA. 

@ USE N. ,NASA. 

@ MAP, IS ,N.VER2 
IN N.PERM 
IN N.FNAL 
IN N.TBAL 
IN N.MPAL 
IN N.CFHE 
IN N.SCZ93 
END 

@ XQT N.VER2 

This deck must be followed by the data deck described in 
Section A-3. In this example the element PERM plus the property elements are mapped 
into an absolute element called VER2. It should be pointed out that once these 
absolute elements have been mapped they are available for subsequent runs 
without being re-mapped. In future runs, if VER2 is called upon to be 
executed, the element consisting of PERM plus is the property elements listed 

above are automatically used. 
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3. Description of Program Units 

This section describes each subprogram which makes up the Permanent 
and Selective Part of the program. 

The description of each subprogram follows the same format. The 
name of the subprogram is given and it is categorized as either a FUNCTI0N 
or SUBR0UTINE subprogram. The argument list is also stated. 

The objective of each of the subprograms is stated. All input 
and output variables in the argument list are defined and units are given 
for all dimensional input and output quantities. 

The calling subprograms and those subprograms called are listed 
for each program unit. 

The program units are subdivided into groups of subprograms each 
of which has a common overall function. For example, all subprograms that are 
used for the calculation of the aerodynamic heating of the radiator panel 
are grouped into a common block of programs. In all, there are 62 program 
units with 42 units on the permanent part and 20 units in the Selective Part. 
The Permanent Part contains 42 program units grouped in seven sets from (a) 
through (g) as follows : 



a) the principal integration program set with 
1. MAIN , the main program unit which 


provides 

°0 input and output 
8) unit conversion 

Y) computation of design parameters 
8) preparations for principal integration, 

calls 

a) RKS , DERIVM, CNTLM, for principal 
integration 

8) FLSTRT to establish the initial flow field 

Y) FUNCTION subprograms for fluid and structural 
material properties ,, . 

<5) QINCID, SHAPEF (see group c) 

e) mathematical procedure subprograms. 

SUBROUTINE DERIVM 

(Y,DY,TIME) which is called from RKS, 

receives 

а) the current system state variables Y 

б) the current time TIME, 

provides 

a) all derivatives of the variables Y with 
respect to time. 

calls 

а) C0NVEC , QRAD 

б) FLSTRT to establish the flow field 

y) FUNCTION subprograms for fluid and 
structural material properties 

6) mathematical procedure subprograms. 

SUBROUTINE CNTLM 

(Y,DY,DX,X,NTRY, IFVD) which 

is called from 

RKS, 

receives 

а) all system state variables Y 

б) their derivatives DY 

✓ 

y) current time X, step size DX, 

provides 

а) output during integration 

б) integration step size control 

Y) integration termination criteria. 

calls 

a) HFL, fluid enthalpy subprogram 
8) mathematical procedure subprograms, 
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b) the secondary integration program set, necessary to establish the 
dynamic fluid flow field, with 

4. SUBROUTINE FLSTRT (RE, PR, DELTA) which 


is called from 


MAIN, DERI VM, 

receives 

a) 

RE, the mean Reynolds number 


6) 

PR, the mean Prandtl number 


Y) 

DELTA, the diameter- to-length ratio for 
the tube. 

provides 

a) 

fluid flow (initial and current) conditions 


6) 

table heading for initial conditions. 

calls 

a) 

TRNSPT to compute friction factor and 
convective film coefficient (in non- 
dimensional form) 


6) 

RKSF, DERIVL, CNTLN 


Y )' 

fluid flow property subprograms 


6) 

mathematical procedure subprograms. 

SUBROUTINE DERIVL 

(Y,DY,X) which 

is called from 


RKSF, 

receives 

a) 

local fluid flow variables Y 


6) 

position along flow channel axis, X, 

provides 

a) 

spatial derivatives DY of Y, 

SUBROUTINE CNTLN 

(Y,DY,DX,X,NTRY, IFVD) which 

is called from 


RKSF, 

receives 

a) 

the flow variables Y 


6) 

their spatial derivatives DY 


Y) 

position X, interval DX along channel axis, 

provides 

a) 

output when initial conditions are 
established 


6) 

termination criteria for RKSF integration. 

SUBROUTINE TRNSPT 

(RE,PRL,DELTA,FR,FNU) which computes non- 


dimensional friction factor and convective film coefficient, 
Eqs. 3.7 through 9 and Eqs. 3.16 through 18, and 

is called from FLSTRT, 

receives «) Reynolds number RE 
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6) 

Prandtl number PRL 


Y) 

diameter-to-length ratio DELTA, 

provides 

a) 

FR = 4f/S, Eq. 3.32 


6) 

FNU = 4 N Nu /(S N^p, Eq. 3.37 


c) The program set used to calculate the incident net radiant flux, 
with, 


8. SUBROUTINE QINCID which . .. 

is called from MAIN, 

receives a) MST0TR control option for simulation of 

steady state (= 1) or transient (= 2) 
conditions 


8) NSRD the number of radiator sides that 1 
are exposed to the environment 

y) PHN the angle between the reference 

used in the MRI file data and the normal 
to the fin panel, 


provides a) data transfer from the MRI Incident 

Radiation Computer program or card input 

8) averaging of radiant flux over circumference 
of tube cross-section, 

9. SUBROUTINE SHAPEF which 

is called from MAIN, 

provides exchange function SS as defined by Eq. 6.15, 

see Appendix D, 

10. SUBROUTINE TRMATX which 


is called from QRAD, 

provides the evaluation of the transfer matrix M. . , 

defined by Eq. 6.7, 

11. SUBROUTINE EXITAV which 

is called from QRAD, 

provides the evaluation of the excitation vector P. 

in Eq. 6.18, 

12. SUBROUTINE ABSpRB which 


is called from 
receives 


QRAD, 

a) total hemispherical emittance e, see Eq. 6.5 

8) temperature distributions on the fin and 
tube surfaces 

y) exchange function SS 


3 
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6 ) 

provides 

calls a) 

e) 

Y) 

6 ) 

13. SUBR0UTINE QRAD which 


functions x. . and x. 
and 6.11, ljk 

adsorptance matrix 


SHAPEF for SS 
EMIT for (T) 


defined by Eqs. 6.10 
as defined by Eq. 6.9, 


AVGEMT for x. . and x, 

ij ijk 

DEFINT for integration, 


is called from 


DERIVM, 

receives 

a) 

transfer matrix and excitation vector 
M„ and , respectively. 

provides 

a) 

inversion of transfer matrix 


8) 

solution of radiosity equations 


Y) 

grid mapping. 

calls 

a) 

TRMATX 


8) 

EXITAV 


Y) 

YINT, INTERP, 


d) the program set used to calculate the aerodynamic heating during 
ascent and reentry, consisting of 


14. FUNCTI0N ALTVEL (TIME, I0PTN , AVA, AVR, TA, TR, NA, NR) which 
interpolates user-supplied altitude and velocity profiles for 
a given instant, and 


is called from 


C0NVEC , 


receives 


а) current real time, TIME, in seconds, 
measured from start of transition phase, 
i.e., ascent: TIME = 0, ALTVEL = AVA (1) 

reentry: TIME = 0, ALTVEL = AVR (1) 

б) ordered pairs of time and altitude 
(TA,AVA) and (TR,AVR) for ascent and 
reentry, respectively 

y) NA and NR, the number of ordered pairs 
and elements in TA,AVA and TR,AVR, 
respectively 
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6) I0PTN = 0 for ascent 
= 1 for reentry, 

provides data interpolation, 

calls YINT for interpolation, 

15. SUBR0UTINE ATM0S (ELEV, TATM, CATM) which computes the atmospheric 
temperature and the speed of sound as functions of altitude and 

is called from C0NVEC, 

receives a) current shuttle elevation, ELEV, in feet 

provides atmospheric temperature, TATM in degrees 

R and speed of sound, CATM, in ft/sec, 

calls YINT for interpolation, P0LY for evaluation 

of polynomials , 

16. SUBROUTINE REFP (ELEV, REFTP, REFPR, REFVIS, REFRH0 , REFK, REFCP, 
REFGAM) which computes the atmospheric air properties as a 
function of elevation and of high speed reference temperature 
and 

is called from C0NVEC, 

receives a) current orbiter elevation, ELEV, in feet 

6) current high speed reference temperature 
REFTP, in degrees R, 

provides atmospheric Prandtl number REFPR 

(dimensionless) , dynamic viscosity REFVIS 
in lbm/ (ft sec), density REFRE0 in lbm/ft 3 , 
thermal conductivity REFK in Btu/(sec ft R) , 
specific heat at constant pressure REFCP 
in Btu/ (lbm/R) and the ratio of specific 
heats REFGAM (dimensionless) , 

calls a) YINT for interpolation 

B) P0LY for polynomial expansion, 

17. SUBROUTINE NUS (MACHN0 , TATM, CATM, TIN, PRATM, VISATM, RH0ATM, 
STAGX, VERTX, NSRAD, NUS1) which provides a single Nusselt number 
for the calculation of the aerodynamic heating of the orbiter 
radiator system at the current time and 


is called from 


C0NVEC , 


receives 

a) 

current orbiter Mach 
(dimensionless) 

number MACHN0 


B) 

current atmospheric 
degrees R 

temperature TATM, in 
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provides 


calls 


y) current atmospheric velocity of sound, 

CATM, in ft/ sec 

e) temperature of the coolant fluid at the 
inlet plane TIN, in degrees R 

C) atmospheric Prandtl number PRATM, 

evaluated at the high speed reference 
temperature (dimensionless) 

n) atmospheric density RH0ATM, evaluated 

at the high speed reference temperature 

in lbm/ ft* 

6) distance from the stagnation point on the 
shuttle vehicle to the midpoint of the 
radiator panel measured along a streamline, 
STAGX, in feet 

O overall dimension of the radiator panel 
measure in a direction parallel to the 
acceleration of gravity, VERTX, in feet 

<) integer value indicating number of non- 
adiabatic sides of the radiator panel, 
NSRAD, i.e., NSRAD = 1 for single non- 
adiabatic surface or NSRAD = 2 for both 
sides of radiator being non- adiabatic, 

mean Nusselt number NUS1, for meteoroid 
layer and fin at current orbiter elevation 
during either ascent or reentry 
(dimensionless) , 

YINT for interpolation, 


18. SUBROUTINE C0NVEC (TIME) which computes normalized aerodynamic 
heating of fin and meteoroid protection layer for both ascent 
and reentry phases of orbiter operation and 


is called from 
receives 


provides 


calls 


DERIVM, 

a) current real time TIME, in seconds, 

measured from start of transition phase; 
i.e., TIME must be measured with respect 
to same reference point as user supplied 
time arrays in FUNCTI0N ALTSH and VELSH, 

normalized convective flux from fin sur- 
face C0NFN (I,J), and from meteoroid 
protection surface C0NMP (I) , 

a) ALTVEL for altitude and velocity of orbiter 

B) ATM0S for atmospheric properties 

Y) REFP for atmospheric properties evaluated 
at high speed reference temperature 
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6) YINT for interpolation 
e) NUS for Nusselt number. 


19. FUNCTION CP AIR (T) which computes the specific heat at constant 
pressure for atmospheric air and 


is called from 

receives 

provides 

calls 


ATM0S, REFTP and C0NVEC , 

absolute temperature T in R, 

specific heat at constant pressure in 
Btu/lbmR, 

P0LY for polynomial expansion. 


20. FUNCTI0N ENTAIR (T) which computes the enthalpy of atmospheric 
air and 


is called from 

receives 

provides 


C0NVEC , 

absolute temperature T in R, 
enthalpy in Btu/lbm, 


21. FUNCTI0N TNH (ENT) which computes the temperature of atmospheric 
air for a given enthalpy and 


is called from 
receives 
provides 
calls 


C0NVEC , 

enthalpy ENT in Btu/lbm, 

absolute temperature in R, 

a) CPAIR for the specific heat at constant 
pressure of atmospheric air 

3) ENTAIR for the enthalpy of air. 


e) the program to compute the thickness of the meteoroid protection 
layer 


22. FUNCTION TK (GAMMA, A, BETA, DENSM, THETA, PHI, AN ALPHA, VELM, 

P0, TAU, DENST, W, TNN, AMAN, TIN, R0UT) which computes the 
thickness of the meteoroid protection layer for given environ- 
mental conditions, tube and manifold areas, experimental constants 
and protection layer properties and 

is called from MAIN, 

receives a) nondimensional experimental constant, 

GAMMA, see Eq. 8.9 
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provides 

calls 


8) nondimens ional experimental constant A, 
see Eq. 8.9 

Y) nondimensional constant which relates 
meteoroid flux to mass, BETA 

<5) density of the meteoroid particle, DENSM, 
in gm/ cm 3 

e) nondimensional experimental constant 
THETA, see Eq. 8.9 

?) nondimensional experimental constant 
PHI, see Eq. 8.9 

n) nondimensional constant AN, used to de- 
scribe penetration depth as a function of 
angle of incidence of meteoroid particle 

9) velocity of meteoroid particle relative 
to radiator panel, VELM, in ft/sec 

O probability of no damage caused by impact 
of meteoroid, P0 

k) time TAU, the radiator panel is exposed 
to meteoroid environment, in days 

A) density of protection material DENST, 
in lbm/ f t 3 

y) axial length of single flow channel exposed 
to meteoroid environment, W, in inches 

v) integer number of flow channels, TNN, 

K) area of the manifold that is exposed, AMAN, 
in ft 2 

4>) temperature of coolant fluid at inlet plane 
of flow channel, TIN, in degrees R 

*0 outside radius R0UT of the unprotected flow 
channel, in inches, 

meteoroid protection thickness TK, in inches, 

ELAS for modulus of elasticity of protection 
material , 


f) the program set which calculates the location of the adiabatic plane 
between two flow channels under conditions of non-symmetrical loading 
of the tubes and for curved radiator panels, consisting of 

23. SUBR0UTINE SHADE which zeros those elements in the array Q which 
correspond to portions of the radiator panel that do not receive 
solar, albedo and planetary irradiation and 



is called from 


TCALC and QINCID, 


receives 


provides 


a) NS, number of flat sides to the polygon 
used as a reference body by MRI program 
for calculating incident radiant flux 

S) radiant flux array Q in Btu/hr ft 2 . 

Array must have NS or greater values, 

zeros in part of the Q array which correspond 
to fin elements in the shade, 


24. SUBROUTINE TCALC which calculates the equivalent sink temperatures 
from incident flux data and. 


is called from MAIN, 

receives NTF the number of flat fin segments. 


provides 

calls 


TSTAR the equivalent sink temperature in R, 
SHADE for incident flux data, 


25. SUBROUTINE EFFICY which evaluates the one- dimensional fin 
efficiency and its derivative as a function of the fin base 
temperature and 


is called from 
receives 


provides 


calls 


TTIPS, 

а) the dimensionless conductance parameter NC 

B) the derivative of the conductance parameter 
with respect to the base temperature NCPR 
in R -1 

y) the ratio of the sink to base temperature 
TS0TB 

б) the fin base temperature TB in R, 

а) the fin efficiency ETA 

б) the derivative of the fin efficiency with' 
respect to the base temperature ETAPR in 

R" 1 , 

P0LY for polynomial expansion. 


26. SUBROUTINE TTIP which computes the tip temperature of a one- 
dimensional fin and 


is called from ADIABH, 

receives a) the dimensionless conductance parameter NC 
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8) the ratio of the sink to base temperature 
TS0TB , 

provides the ratio of the difference between the fin 

tip temperature and sink temperature to the 
difference between the base and sink temper- 
ature, 

calls P0LY for polynomial expansion, 

27. FUNCTI0N NUSA which evaluates the Nusselt number for the coolant 
fluid flow for both laminar and turbulent flow and 


is called from 


ADIABH and TTIPS, 


receives 

a) 

the Reynold number for the 

fluid REY 


8) 

the Prandtl number for the 

fluid PR 


Y) 

the ratio of the internal 
tube length D0L, 

tube diameter to 

provides 


the fluid Nusselt number. 



28. SUBR0UTINE TTIPS which uses the Newton-Raphson method to numerically 

solve for the base temperature of the one-dimensional fin and 

is called from ADIABH, 

receives a) the integer IT used to specify the tube 

number for which the fin tip temperatures 
are calculated. Numbering system is 
consistent with the input parameters TIN 
(inlet temperature), PIN (inlet pressure), 

T (fin thickness), AL (tube length) and M 
(mass flow rate) 

8) the integer ITST which designates the sink 
temperature data stored in the array TSTAR 

y) the distance Hi to the adiabatic plane for 
the fin attached to the left of the tube 
number IT, in inches 

<5) the distance H2 to the adiabatic plane for 
the fin attached to the right of the tube 
number IT, in inches, 

provides a) the adiabatic plane temperature for the 

fin attached to the right of tube IT in R 

8) the adiabatic plane temperature for the 
fin attached to the left of tube IT in R, 

calls a) NUS for the coolant fluid Nusselt number 

8) EFFICY for fin efficiency 

Y) TTIP for the fin tip temperature. 



29. SUBR0UTINE ADIABH which determines the position of the adiabatic 
planes so that there is a continuity in adiabatic plane tempera- 
ture of the one- dimensional fin and so that the heat convected 
from the fluid equals the heat radiated from the fin and 


is called from MAIN, 


receives 


provides 


calls 


the integer ITST which designates sink 
temperature data stored in the array TSTAR,. 

printout of table headings and system 
parameters when radiator is non- symmetric ally 
loaded and diagnostic printout if distance to 
adiabatic plane exceeds the spacing between 
the tubes. For this case one tube is trans- 
ferring heat to the fin while the adjacent 
tube is gaining heat. TTIPS also provides 
the distances to all the adiabatic planes 
within the non-symmetrical radiator section, 

а) property subroutines for fin and fluid 
properties 

б) NUS for the coolant fluid Nusselt number 

y) TTIPS for the tip temperature of the fin 

6) PDERIV for determination of the elements 

of the matrix P. . 

ij 

e) FMINV for matrix inversion. 


30. FUNCTION PDERIV which calculates the partial derivative 36^/ 9h^ 
(see Section 11-14 of the System Analysis Manual) and 


is called from 
receives 


provides 

calls 


ADIABH, 

a) the integer J denotes the tube under 

consideration, J = -1 denotes tube IT - 1; 

J = 0 denotes tube IT and J = +1 denotes 
tube IT + 1 

8) the integer IT which specifies the tube 
number to be considered 

y) the integer ITST which designates sink 

temperature data stored in the array TSTAR, 

the partial derivative 3h^/3hj 

TTIPS for evaluation of the temperature at 
the adiabatic plane. 
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g) the program set consisting of all the mathematical procedures 
discussed in detail in Part III of the System Analysis Manual' 

31. SUBROUTINE RKSF (see Section III. 15) which 

is called from FLSTRT, for dynamic fluid flow field, 

calls DERIVL , CNTLN , 

32. SUBROUTINE RKS (see Section III. 15) which 

is called from MAIN, for principal integration, 

calls DERIVM, CNTLM, 

33. FUNCTION PpLY (X,A,M) (see Section III. 16), 

34. FUNCTION YINT (X,Y,M,N,P) (see Section III. 17), 

35. SUBROUTINE DDX (Y,DY,DX,N) , 

(see Section III. 18) 

36. SUBROUTINE D2DX2 (Y ,D2Y,DX,N) , 


37. FUNCTION DEFINT (Y,DX,N), 

38. SUBROUTINE FINT (Y,YO,DX,N,F) , 


(see Section III. 19) 


39. SUBROUTINE INTERP (NX1 ,MZ1 ,XX1 , ZZ1 ,YY1 ,NX2 ,MZ2 ,XX2 , ZZ2 , YY2) maps 
a two-dimensional function YY1 from one grid (XXI, ZZ1) onto 
another grid (XX2,ZZ2) and 


is called from QRAD, 

receives a) the number of nodal points (NX1,MZ1) of 

the original grid 

B) the function YY1 

y) the number of nodal points (NX2,MZ2) of 
new grid, 

provides the function YY2 at the new grid, 

calls YINT, 

is restricted to NX1 ,MZ1 ,NX2 ,MZ2 <. 10, 


40. FUNCTION DEFNT serves to integrate by the trapezoidal rule and 
is used where random error accumulation is more critical than 
truncation errors, and V accepts through its argument list 
(Y, DX, N) 



a) the array Y of N equally spaced ordinates 

P) the interval DX 

y) the integer constant N, representing the number of nodal 
points. 

returns the integral of YCS) as DEFNT. 

41. SUBROUTINE MTXINV (N,M) performs two related tasks: 

(i) to solve a system of N linear algebraic equations, 
when M = N + 1, 

(ii) to invert an N x N invertible matrix, when 
M = 2N. 

For task (i) it 

accepts through C0MM0N TRMTX (k,l) the augmented coefficient 

matrix (see Chapter III. 7). Here k. and 1 represent 
integer constants, k=N, 1 = N + 1. 
through the argument list (N,M) the rank N of the 
coefficient matrix and M = N + 1, two integer 
constants. 

returns the solution in TRMTX (I,J), 1=1, ..., k, J = k. 

For task (i) it 

accepts through C0MM0N TRMTX (k,l) 

the k x k matrix in the first k columns of 
TRMTX (I,J), 1= 1, ..., k; J= 1, ..., k and the 
k x k identity matrix in the second k columns of 

TRMTX (I,J) 1 = 1 k; J = k + 1, . . . , 2k. 

Hence 1 = 2k. 

through the argument list (N,M) the rank N of the 

coefficient matrix and M = 2N, both integer constants 

returns the inverted matrix in the second k columns of TRMTX 
(I,J), 1=1, ..., k; J = k + 1, ..., 2k. 

42. SUBR0UTINE FMINV serves to solve N linear algebraic equations 

by the method used for SUBR0UTINE MTXINV. It's calling provisions 

are designed to accommodate the requirements of ADIABH. It 
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receives, through its argument list (A,X,N,M) 

o') the (N x N) , invertible coefficient matrix A 
3) the know, N-dimensional vector X 
y) the rank (N) of A, an integer constant 
6) M = N + 1 

where N ^ 25 

returns the solution Y = A"^X placed in the array X 
is called from ADIABH. 
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The Selective Part contains five groups of thermophysical properties, with 
the total of 20 program units. Each group constitutes a property package, 
one each for 

(i) the coolant fluid (8 program units) 

(ii) the flow channel material (3 program units) 

(iii) the meteoroid protection layer (4 program units) 

(iv) the fin material (3 program units) 

(v) the thermal control coating (2 program units). 

The choice of a particular material combination must be reflected in the 
corresponding data specification of the Input Data Set, (see Section 3f of,, 
Chapter A) • 

a) The coolant property subprogram set consists of eight elements, 43 

through 50 

43. FUNCTI0N RH0F (P,T) which computes the fluid density RH0F in 
slug/ft 3 as a function of pressure P in lbf/ft 2 , and absolute 
temperature T in R, 

44. FUNCTI0N BETA (RH0,T) which computes the isobaric expansion 
coefficient 8, Section II.C.ll, in 1/R as a function of density 
RH0 in slug/ft 3 , and temperature T in R, 

45. FUNCTI0N CAPPA (P,T) which computes the isothermal compressibility, 
Eq. 11.3, in ft 2 /lbf, as a function of pressure P in lbf/ft 2 , and 
temperature T in R, 

46. FUNCTI0N CPF (RH0,T) which computes the specific heat at constant 

pressure CPF, in Btu/(slug R) , as a function of density RH0 in 
slug/ft 3 , and temperature T, in R, . , 

47. FUNCTI0N HFL (RH0,T) which computes the fluid enthalpy HFL, in 
Btu/slug, as a function of density RH0 in slug/ft 3 , and tempera- 
ture T, in R, 

48. FUNCTI0N VISC (RH0,T) which computes the fluid dynamic viscosity 
VISC, in slug/ (ft sec) as a function of density RH0 in slug/ft 3 , 
and temperature T, in R, 
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49. FUNCTI0N THCF (RH0,T) which computes the fluid thermal conductivity 
THCF in Btu/ (hr ft R) as a function of density RH0 in siug/ft 3 , 
and temperature T, in R, 

50. FUNCTI0N PF (RH0,T) which computes the fluid pressure PF in lbf/ft 2 
as a function of density RH0 in slugs/ft 3 , and absolute temperature 
T, in R. For specific details of the coolant fluid properties, 
see Section II. 11 and 12 and Appendix B of the System Analysis 
Manual. 

b) The coolant channel material property subprogram set consists of 

three elements, 51 through 53. 

51. FUNCTI0N THCTB (T) computes the thermal conductivity THCTB in 
Btu/ (hr ft R) , of the tube wall material as a function of the 
absolute temperature T, in R, 

52. FUNCTI0N DTHCTB (T) computes the relative change of. thermal 
conductivity k with temperature T, namely (dk w /dT)/k w in 1/R, 
as a function of temperature T, in R, 

53. FUNCTI0N CPTB (T) computes the specific heat at constant pressure 
CPTB, in Btu/ (slug R) as a function of T, in R. 

For specific details of the coolant channel properties, see Appendix A of the 
System Analysis Manual. 

c) the meteoroid protection layer material properties are coded in a set 

of four elements, 54 through 57 

54. FUNCTI0N THCMP (T) , 

55. FUNCTI0N DTHCMP (T) , 

56. FUNCTI0N CPMP (T) , 

57. FUNCTI0N ELAS (T) computes the modulus of elasticity ELAS in 
lbf/in 2 as a function of temperature T, in R. 

For specific details of the meteoroid protection properties, see Appendix A 
of the System Analysis Manual. 

d) The fin material properties required are coded in three FUNCTI0N sub- 
programs, numbered 53 through 60 
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58. FUNCTI0N THCFN (T) computes the thermal conductivity THCFN, in 
Btu/(hr ft R) for the fin material as a function of absolute 
temperature T, in R, 

59. FUNCTI0N DTHCFN (T) computes the relative change of thermal con- 
ductivity, namely (dk^/dT)/k^ in 1/R, as a function of temperature 
T, in R, 

60. FUNCTI0N CPFN (T) computes the specific heat at constant pressure 
CPFN, in Btu/(slug R) as a function of temperature T, in R. 


For specific details of the fin material properties, see Appendix A of the 
System Analysis Manual. ' 


e) Optical properties of the thermal control coating are coded in two 
subprograms, numbered 61 and 62 

61. FUNCTI0N EMIT (T) computes the total hemispherical emittance in 

accordance with Eq. 6.5 of the System Analysis Manual as a function 
of the surface temperature T in R, and is called by ABS0RB and 
EXITAV, 


62. SUBR0UTINE AVGEMT (T,XX,XXX,N) evaluates the expressions for x^ 
and given as Eqs . 6.10 and 11 in the System Analysis Manual 

for N interacting surface elements. It is called by ABS0RB (see 
program No . 12) . 


For specific details of surface coating properties, see Appendix C of the 
System Analysis Manual. 

This completes the Source Deck discussion. Recall that the Source 
Deck composition implies a particular selection of coolant fluid and structural 
materials . 





117 


4. Program Listing 

This section contains the entire program listing in alphabetical 
order of subprogram name. 

Labelled common blocks and external program references and 
printed out prior to the listing of each program unit. 

All program units have been documented with comment cards for 
the convenience of the user. 

All program units are summarized in a table of contents which 
follows the final unit listing. The table of contents lists all program 
elements as either symbolic or relocatable and it also prints the time and 
date when the program element was most recently stored on FASTRAND drum. 
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